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INTRODUCTION 


Combustion  of  the  liquid  fuels,  through  atomization  and  vaporization,  will 
remain  the  primary  mode  of  fossil  fuel  energy  conversion  for  the  next  several 
decades.  It  finds  a  wide  range  of  applications  such  as  in  variety  of  gas  turbine 
engines,  diesel  and  fuel-injected  spark  engines,  industrial  boilers  and  liquid  rocket 
motors.  The  spray  combustion  in  all  these  systems  is  a  complex  phenomena 
involving:  (i)  atomization  of  liquid  fuels,  which  produces  droplets  of  different  sizes 
and  velocities  (ii)  mixing  of  the  droplets  with  oxidant  and  hot  combustion  products 
(iii)  vaporization  of  the  droplets  (iv)  mixing  of  the  fuel  vapor  with  surrounding  gas  by 
diffusion  and  convection  processes  and  (v)  chemical  reaction.  The  tasks  of 
improving  combustion  efficiency  and  pollutant  levels  require  a  fundamental 
understanding  of  these  physical  and  chemical  processes.  Because  of  the  complexity 
of  the  problem,  a  simultaneous  enhancement  is  needed  in  experimental  measurements 
and  predictive  capabilities  of  spray  combustion  phenomena.  With  this  objective  in 
mind,  an  ARO-sponsored  theoretical-experimental  program  has  been  underway  at 
Carnegie-Mellon  University  for  the  past  four  years. 
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THEORETICAL  EFFORT 

This  section  summarizes  the  theoretical  and  numerical  efforts  in  enhancing  the 
fundamental  understanding  of  spray  combustion  processes.  This  research  has  been 
supported  by  Army  Research  Office  and  covers  a  four-year  period  from  1979  to 

1983.  During  this  period,  several  major  tasks  have  been  completed.  These  include: 

1.  the  study  of  enclosed,  unsteady  one-dimensional  laminar  and  turbulent 
flame  propagation  in  homogeneous  and  heterogeneous  mixtures, 

2.  the  efficiency  comparison  of  several  explicit  and  implicit  finite-difference 
schemes  for  solving  laminar  homogeneous  and  heterogeneous  enclosed 
deflagrations, 

3.  the  efficiency  and  accuracy  comparison  of  explicit,  implicit  and  operator¬ 
splitting  techniques  for  axisymmetric  spray  vaporization  problem, 

4.  the  study  of  one-dimensional  steady  laminar  flame  propagation  in  one- 
phase  and  two-phase  mixtures, 

5.  the  development  of  hybrid  (Eulerian-Lagrangian)  numerical  code  for  one¬ 
dimensional,  reacting  air-fuel  droplet-fuel  vapor  flows, 

6.  computation  of  spray  flame  speeds  and  a  preliminary  comparison  with 
experimental  results, 

7.  a  parametric  study  to  examine  the  effects  of  droplet  size,  equivalence 
ratio,  fuel  type,  droplet  heating  models,  chemical  kinetics,  and  grid  size 
and 

8.  the  examination  of  different  gas  and  liquid  phase  models  which  are  used 
in  spray  vaporization  and  combustion  processes. 

Task  1  is  discussed  in  a  paper  (see  addenda)  presented  at  the  1980  ASME 
Winter  Annual  Meeting.  In  this  study,  the  gas  turbulence  was  modeled  through  a 
two-equation  (k-t)  model  and  the  particle  turbulence  was  neglected.  The  major 
conclusion  was  that  the  characteristics  of  turbulent  flame  propagation  in  both 
homogeneous  and  heterogeneous  mixtures  depend  on  a  turbulence  delay  time  defined 
by  the  parameter  k/c,  and  turbulence  generation  time  which  is  related  to  the  chemical 
reaction  time.  Task  2  evaluated  the  effectiveness  of  several  explicit  and  implicit 
finite-difference  methods  for  computing  the  unsteady  laminar  flames  in  homogeneous 
as  well  as  heterogeneous  mixtures.  The  details  are  given  in  Refs.  2  and  3  (see 
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addenda).  A  similar  study  has  been  performed  under  Task  3,  where  vaporizing  spray 
is  considered  in  a  gaseous  jet  flow.  The  effectiveness  curves  are  obtained  for  four 
different  finite-difference  schemes.  The  results  are  detailed  in  the  paper  submitted 
for  publication  in  Journal  of  Computational  Phvsics  (Ref.  4). 

These  are  two  major  objectives  for  Tasks  4-6.  The  first  is  to  develop  an 
efficient  numerical  code  to  model  the  two-phase  processes  which  occur  in  laminar 
spray  flames.  The  results  of  this  code  are  expected  to  provide  global  as  well  as 
detailed  local  information  about  the  spray  flame  phenomena.  The  global  information 
may  include  the  flame  speed,  the  influence  of  various  parameters  on  the  flame 
speed,  and  the  identification  of  dominant  burning  mechanism  of  the  spray  flame. 
The  second  objective  is  to  provide  numerical  results  for  making  comparison  with  the 
results  of  the  parallel  experimental  study.  A  hybrid  Eulerian-Lagrangian  code  has 
been  developed,  the  laminar  spray  flame  speeds  have  been  computed,  a  limited 
parametric  study  has  been  conducted,  and  some  preliminary  comparison  has  been 
made  with  the  experiments.  These  efforts  are  detailed  in  Ref.  5-7.  At  present  the 
numerical  code  is  being  modified  to  include  the  polydispersed  spray,  more 
sophisticated  droplet-heating  models  and  the  multicomponent  fuel  droplets.  The 
probabilistic  character  of  the  spray  and  the  droplet  group  combustion  behavior  will 
also  be  included  in  a  later  study. 

Under  Task  7,  the  hybrid  numerical  code  has  been  employed  to  perform  a 
detailed  parametric  study  of  flame  propagation  in  confined  sprays.  Some  results  are 
given  in  Ref.  8.  Additional  results  will  are  presented  in  Ref.  9.  Under  Task  9, 
capabilities  have  been  developed  to  include  more  sophisticated  gas-phase  and  liquid- 
phase  models  for  representing  the  inter-phase  processes.  A  comparison  of  various 
models  for  an  isolated  droplet  and  for  a  vaporizing  spray  is  given  in  Ref.  10. 
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EXPERIMENTAL  EFFORT 

An  experimental  study  of  ignition  and  flame  propagation  in  sprays  is  achieved 
by  studying  evaporation  and  combustion  of  a  well  controlled  spray.  During  this 
research  a  unique  spray  generator  is  designed  and  developed  which  is  capable  of 
producing  a  well  controlled  droplet  flow.  Later  this  spray  is  burned  and  effect  of 
important  parameters  such  as  droplet  diameter,  air/fuel  rates,  and  fuel  type  on  the 
combustion  phenomena  is  studied.  Also,  an  analytical  study  of  the  radiative 
evaporation  of  the  fuel  droplet  flow  at  the  upstream  of  the  flame  is  done.  In 
following  sections  a  brief  discussion  of  each  of  above  work,  with  reference  to 
corresponding  publications  are  given. 

SPRAY  GENERATOR 

In  order  to  generate  a  well  controlled  droplet  flow  spray,  a  multi-orifice 

impulsed  spray  generator  with  two  different  modes  of  operation  is  developed.  The 
impulse  jet  operation  mode  covers  a  wide  range  of  droplet  frequencies,  and  produces 
arrays  of  droplets  with  controllable  spacing,  but  with  a  wide  size  spectrum.  This 
ejector  also  operates  in  pressurized  flow  mode.  In  this  mode,  optimum  disturbance 
frequency  breaks  up  the  liquid  jets  into  uniform  size  droplets,  resulting  in  a 
monosized  spray.  Precise  control  of  fuel  flow  rate  and  droplet  size  is  obtained. 
The  details  of  design  and  operation  of  multi-orifice  impulsed  spray  generator  is  given 
in  publications  generated  from  this  research  listed  in  references  11-13.  Reference  13 
is  included  in  the  appendix. 

SPRAY  COMBUSTION  STUDIES 

In  order  to  study  the  combustion  behavior  of  the  spray  generated  by  the 

present  device,  a  special  experimental  set  up  is  fabricated  (13).  The  experimental 

study  of  the  monodispersed  spray  of  various  fuels  identified  three  major  combustion 
phenomena.  The  all-blue  flame  of  premixed  gas  combustion;  the  isolated  droplet 

burning  and  propagation;  and  the  mixed  feature  of  premixed  gas  combustion  and 
droplet  burning  where  droplet  cluster  combustion  is  in  this  category. 
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Specific  parametric  studies  have  revealed  that  a  dominant  factor  which  is 
affecting  the  spray  combustion  phenomena  and  mechanism  is  the  vapor  fuel  to  liquid 
fuel  ratio  at  the  flame  front.  The  amount  of  vapor  fuel  is  dependent  upon  the 
droplet  vaporization  process  which  is,  in  turn  dependent  upon  order  factors.  The 
mechanism  of  the  flame  propagation  in  sprays  is  governed  by  two  different 
phenomena.  One  is  kinetically  controlled  phenomena,  and  the  other  is  diffusion 
controlled  one.  If  there  is  enough  prevaporization,  or  there  is  enough  fuel  vapor  in 
the  fresh  mixture  so  that  the  fuel  vapor/air  mixture  ratio  ahead  of  the  ignition  zone 
is  flammable  then  flame  propagation  is  governed  by  the  premixed,  laminar  or 
turbulent  gas  flame  propagation.  But  if  there  is  not  enough  vaporization,  then  overall 
fuel  vapor/air  mixture  would  not  be  flammable,  although,  it  is  observed  that  flame 
propagates  in  such  sprays  also.  This  type  of  flame  propagation  is  different  than  that 
of  premixed  gas  combustion.  In  such  conditions  the  spray  is  either  very  dilute  that 
relay  type  flame  propagation  is  observed,  where  flame  of  an  ignited  droplet 
increases  in  size  due  to  oxygen  retardation,  and  eventually  ignites  the  adjacent 
droplets.  And,  if  spray  is  very  rich,  the  group  combustion  phenomena  is  observed, 
where  its  mechanism  is  yet  to  be  understood  better.  Both  of  these  phenomena  are 
diffusion  controlled. 

In  order  to  see  the  effect  of  fuel  vapor  fraction  on  the  flame  propagation 
mechanism,  experiments  are  done  with  three  different  fuels.  Hexane,  Heptane,  and 
Octane.  A  multi-orifice  plate  with  47  orifices  with  150/im  diameter  were  used  to  do 
the  experiment.  The  average  initial  droplet  size  for  all  three  cases  was  320  y/m  and 
overall  equivalence  ratio  was  kept  constant  at  1.8.  Applied  frequency  was  3400  Hz. 
Initial  fuel  temperature  was  varied  to  change  the  fuel  vapor  fraction  in  the  spray. 
As  the  initial  fuel  temperature  is  increased,  more  fuel  is  vaporized.  As  a  result  the 
fuel  vapor  fraction  at  the  ignition  point  increases  as  fuel  initial  temperature 
increases.  Figure  _1  shows  a  plot  of  flame  speed  vs.  initial  fuel  temperature.  It  is 
seen  that  the  flame  speed  increases  as  the  fuel  temperature  and  thus  the  fuel  vapor 
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fraction  increases.  But  the  rate  of  this  increase  is  dependent  on  type  of  fuel.  As 
the  fuel  becomes  less  volatile  the  rate  of  increase  of  flame  speed  becomes  less. 
The  flame  speed  for  volatile  fuel  of  iso-hexane  can  go  as  low  as  20  cm/ sec  and 
increase  to  80  cm/sec  for  initial  temperature  of  43°C.  The  results  of  this  work  will 
be  published  in  near  future. 

RADIATIVE  HEATING  OF  THE  DROPLETS 

The  heating  and  evaporation  of  the  fuel  droplet  flow  at  the  upstream  of  the 
flame  front  are  analyzed  considering  the  effect  of  thermal  radiation  from  the  hot 
zone  of  the  combustion  system.  The  radiative  transport  in  a  typical  cylindrical 
combustion  chamber  is  calculated  considering  the  hot  wall,  the  emitting  and 
absorbing  combustion  products,  and  the  absorption  of  the  two  phase  droplet  flow. 
Analytical  solution  is  obtained.  The  heating  and  evaporation  of  a  single  droplet  by 
thermal  radiation  in  cold  air  has  also  been  studied.  The  transfer  number  B  is 
modified  to  include  the  effect  of  radiation.  With  the  presence  of  radiation  the 
change  of  the  droplet  size  will  not  follow  the  simple  d2  rule.  At  the  extreme  case 
of  quiescent  environment  the  droplet  surface  temperature  passes  through  a  peak 
before  the  droplet  evaporates  completely.  As  an  approximation  the  radiative  flux  at 
the  upstream  of  the  flame  front  may  be  considered  as  a  constant,  and  the 
subsequent  calculation  will  be  easier.  Following  the  methodology  outlined  in  this 
work,  complicated  situations  may  be  analyzed.  Details  of  this  work  are  given  in 
publications  generated  under  this  grant  which  are  listed  in  references  (14-15)  and  are 
included  in  the  appendix. 
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One-Dimensional  Turbulent  Flame 
Propagation  in  an  Air-Fuel  Droplet 
Mixture 

The  turbulent  flame  initiation  and  propagation  through  an  air -fuel  drop  fuel  spray 
system  is  numerically  modeled  The  present  work  forms  an  extension  of  the  laminar  flame 
propagation  model  ( 1 ),  where  a  dilute  spray  of  a  single  component  fuel  in  a 
one -dimensional  closed  combustor  is  considered.  The  model  employs  an  unsteady 
formulation  of  the  problem  and  allows  unsteady  mass ,  momentum  and  heat  transfer 
between  the  phases.  The  effect  of  turbulence  on  the  gas  phase  is  incorporated  by  using  the 
time -averaged  governing  equations  with  the  correlation  terms  being  modeled  through  a 
gradient  transport  hypothesis.  The  time -dependent  profiles  of gas  and  liquid  properties  as 
well  as  of  the  turbulent  quantities  are  htained  for  homogeneous  (pre-mixed)  and 
heterogeneous  mixtures.  The  results  indicate  that  the  characteristics  of  flame 
propagation  in  both  homogeneous  and  heterogeneous  mixtures  are  sensitive  to  the  initial 
value  of  the  k/t ,  which  defines  a  turbulence  decay  time ,  and  to  the  chemical  kinetics , 
which  defines  a  turbulence  generation  time. 


NOMENCLATURE 
Cu  Constant 

Cp  Average  specific  heat  at  constant  pressure 

Dt  Eddy  diffusivity 

D  Laminar  mass  diffusivity 

Dp  Particle  diffusivity  n  w 

K  Variable  defined  as  Tpl  '  J 

k  Turbulent  kinetic  energy 

L*  Total  heat  transfer  per  unit  volume  to  a  droplet 
.  of  group  i 

nty  Mass  vaporization  rate  from  droplet  belonging  to 
group  i 

n  Number  density  of  droplets  belonging  to  group  i 
0  Order  of  magnitude 

p  Pressure 

Q  Heat  released  per  unit  mass  of  fuel  burnt 
R  Radius  of  droplets 

T  Temperature 

tp  Particle  response  time 

tt  Turbulent  time  scale 

t  Time  variable 

Ts  Droplet  surface  temperature 

u  Time  averaged  gas  phase  velocity 

ft  Instantaneous  gas  phase  velocity 

Droplet  velocity 

Time  rate  of  change  of  the  fuel  mass  due  to 
chemical  reaction 
x  Space  variable 

Mass  fraction  of  species  i 


GREEK  SYMBOLS 

P  Defined  as  (1  -  Y)/Y 

Y  Specific  heat  ratio 

t  Turbulent  kinetic  energy  dissipation  rate 
Ax  Grid  spacing  in  the  axial  direction 

At  Grid  spacing  in  the  time  direction 

ut  Turbulent  viscosity 

u  Laminar  viscosity 

c k  Constant 

o£  Constant 

\>  Stoichiometric  ratio  of  fuel  mass  to  oxidizer 
mass 

o  Density 

SUBSCRIPTS 

f  Referring  to  the  fuel 

i,j  at  i*h  grid  point  in  x-direction  and  j**1  time 
step 

l  Referring  to  the  liquid  phase 

N  Referring  to  the  neutral  species 

0  Referring  to  the  oxidizer 

S  Droplet  surface 

SUPERSCRIPTS 

1  Fluctuating  part  of  a  variable 

^  Instantaneous  value  of  a  variable 
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The  theoretical  study  of  transient,  turbulent, 
two-phase  flow  in  real  situations  often  leads  to  a 
complicated  set  of  coupled  nonlinear  partial 
differential  equations.  A  practical  solution  to 
these  equations  usually  involves  some  simplifying 
assumptions  as  well  as  the  need  for  numerical  modeling. 
This  approach  has  been  followed  by  several  investi¬ 
gators  ( -C  *  3 , 4 )  Following  the  same  approach,  a  model 
of  unsteady  laminar  flame  initiation  and  propagation 
through  an  air-fuel  drop-fuel  spray  in  a  one- 


dimensional  closed  combustor  was  recently  presented 
A  dilute  poly-disperse  spray  of  a  single  component 
fuel  in  a  mist  form  was  considered.  The  interaction 
between  the  droplets  was  neglected.  The  other 
important  features  included  were:  transient  heating 
and  prevaporization  of  the  droplets,  lag  between  gas 
and  liquid  velocities,  the  effect  of  convection  by 
semi -empirical  relations,  the  coupling  between  droplet 
vaporization  and  environmental  conditions,  and  the 
use  of  a  droplet  ignition  criterion  based  on  the 
concept  of  ignition  Damkohler  number  as  derived  by 
Law(S).  The  initial  gas  and  liquid  phase  properties 
were  specified  in  an  ad  hoc  but  consistent  manner. 

The  ignition  was  simulated  by  an  imposed  increase  of 
temperature  and  an  associated  decrease  in  fuel  and 
oxidizer  concentrations  in  a  small  zone  at  one  end 
of  the  combustor.  Under  appropriate  conditions,  a 
flame  was  found  to  establish  and  to  propagate  across 
the  combustor.  All  the  details  and  the  results  of 
the  laminar  flame  propagation  model  can  be  found  in 
reference  (1). 

The  laminar  flame  model  is  now  extended  to  in¬ 
clude  the  effect  of  turbulence.  The  governing  equa¬ 
tions  used  for  the  laminar  case  are  time-averaged 
which  introduces  the  undetermined  correlation  terms. 

If  the  density  fluctuations  are  disregarded,  these 
correlations  are  of  second  order.  In  order  to  solve 
this  closure  problem,  a  gradient  transport  hypothesis, 
as  suggested  by  Boussinesq ,  is  used.  The  resulting 
turbulent  diffusivity  is  determined  by  employing  the 
two-equation  model  of  turbulence (6) ,  i.e.,  from  the 
solution  of  transport  equations  for  the  turbulent 
kinetic  energy  (k)  and  the  energy  dissipation  rate 
(e) .  A  similar  treatment  of  the  instantaneous  liquid 
phase  equations  will  introduce  the  concept  of  droplet 
or  particle  diffusion.  It  is  hypothesized  that  the 
ratio  of  gas  phase  diffusivity  and  particle  diffusi¬ 
vity  (which  is  like  a  Prandtl  number  for  the  particle 
diffusivity)  depends  solely  on  the  ratio  of  a  particle 
response  time  and  a  characteristic  turbulent  eddy 
time.  Using  this  hypothesis  and  the  Stoke' s  drag  law, 
the  particle  diffusivity  is  found  to  be  a  function  of 

(a)  particle  mass  density/gas  phase  mass  density, 

(b)  square  of  the  ratio  of  the  particle  size  and  the 
turbulent  length  scale  and  (c)  turbulent  diffusivity/ 
laminar  diffusivity.  It  is  argued  that  if  the  turbu¬ 
lent  diffusivity  is  two  or  more  orders  of  magnitude 
higher  than  the  laminar  diffusivity,  the  particle 
diffusivity  may  be  neglected.  For  the  present  calcu¬ 
lations  also,  the  particle  turbulence  is  neglected 
though  its  inclusion  is  straightforward. 

As  the  equation  for  k  would  indicate,  combustion 
is  the  only  turbulence-generating  mechanism  in  the 
present  model.  The  influence  of  turbulence  on  vapo¬ 
rization  and  combustion  is  considered  only  in  an 
indirect  way;  i.e.,  the  turbulence  changes  the 
transport  rates  which  will  influence  the  local  values 
of  the  mean  flow  variables  and  thus  will  influence 
the  vaporization  and  chemical  reaction  rates.  To 
include  the  direct  effect  of  turbulence,  one  needs 
to  consider  how  the  chemical  reaction  and  vaporiza¬ 
tion  rates  are  influenced  by  the  turbulence-induced 
fluctuations  in  composition,  temperature,  and  density. 

One  could  indicate  several  discrepancies  in  the 
use  of  this  type  of  modeling  in  flows  involving 
turbulent  combustion.  For  example,  the  assumptions 
of  gradient  transport,  of  incompressible  turbulence 
behavior  and  of  unimportance  of  direct  effects  of 
turbulence  on  reaction  and  vaporization  rates  remain 
open  to  question (?»8) .  The  major  justification  for 
our  approach  is  that  the  problem  at  hand  is  made 
amenable  to  a  practical  treatment.  Moreover,  the  use 


of  k  -  e  model  has  been  used,  rather  successfully,  in 
the  calculation  of  single-phase  combusting  turbulent 
flows  in  complex  geometries (9 , 10 , 1 1) . 

The  derivations  of  the  governing  equations,  the 
solution  procedure  and  the  numerical  results  are  dis¬ 
cussed  in  the  following  sections.  Since  the  present 
work  is  an  extension  of  the  laminar  flame  propagation 
model  (H,  reference  (1)  has  been  very  frequently  cited 
for  detailed  documentation. 


THE  GOVERNING  EQUATIONS 


The  assumptions  used  in  reference  (1)  are  also 
applicable  to  the  present  case.  If  one  considers  the 
instantaneous  flow  variables,  one  can  use  the  same 
set  of  gas  phase  and  liquid  phase  equations  as  given 
in  reference  (1).  However,  the  instantaneous  flow 
variables  are  now  decomposed  into  a  time-dependent 
mean  and  a  fluctuating  quantity  as: 
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The  assumption  of  incompressible  turbulent 
behavior  has  been  used  in  the  above  equations.  These 
equations  are  substituted  into  the  governing  equations 
of  reference  (1)  and  a  long-time  average  of  the 
resulting  equations  is  taken.  As  a  result  several 
correlation  terms  appear  in  the  equations.  The 
velocity  correlations,  with  the  mean  flow  being 
assumed  one-dimensional,  are  modeled  as 
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The  summation  of  (2),  (3),  and  (4)  gives 
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The  correlations  involving  the  fluctuations  of  a 

scalar  variable  are  modeled  by  using  a  gradient  & 

hypothesis,  i.e.  , 
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where  is  a  turbulent  Prandt 1/Schmidt  number  and  is 
assumed  to  be  unity.  The  eddy  viscosity  ut  is  given 
by 

Ut  -  Cuok2/e  (?) 


> 


2 


I 


1 


Using  the  above  relations,  the  equations  for  Y^,  Yq, 
Yr  and  K  may  be  generalized  into  a  single  equationas 


1±  ♦  u  it  -  i  2_  (PD  »t) 

3t  3x  p  3X  w  •  3x' 


- 1L  S*  (8) 

0  0 


where  $  may  represent  any  of  the  aforementioned  vari¬ 
ables.  It  may  be  mentioned  that  T  has  been  replaced 
by  a  variable  K.  This  eliminates  the  pressure  term 
in  the  energy  equation  (reference  l) .  S$  represents 
the  source  term  and  is  given  in  Table  1  for  each  of 
the  variables. 
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is  the  effective  diffusivity,  given  by 
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It  is  important  to, note  that  a  term  due  to 
droplet  vaporization  (nm  /p)  appears  as  a  sink  term 
in  k  and  c  equations.  This  is  because  the  equations 
are  written  in  a  non- conservation  form.  This  term 
will  be  absent  in  the  conservation  form  of  k  and  e 
equations,  which  means  that  there  is  no  turbulence 
associated  with  the  vaporization  process.  However, 
along  a  streamline  k  and  e  will  change  due  to  vapori¬ 
zation. 

It  may  also  be  mentioned  that  the  term  -2/3  k  • 
3u/3x  appearing  in  these  equations  represents  the 


effect  of  compressibility.  For  example,  in  case  of 
expansion  >  0,  and  the  effect  of  this  term  is  to 
reduce  the^rate  of  increase  of  turbulent  kinetic 


energy.  The  values  of  the  constants  appearing  in  the 
above  equations  are  taken  to  be (19) 
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PARTICLE  DIFFUSION 


Gas  phase  turbulence  can  cause  droplets  or 
particles  to  fluctuate  about  some  mean  velocity. 

This  gives  rise  to  the  concept  of  particle  diffusion, 
which  may  be  approximated  by  adding  a  diffusion  term 
in  the  liquid  phase  equations.  This  term  may  be 
modeled  as 
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The  expressions  for  nty,  w^f  and  L'  are  assumed 
to  remain  unchanged  for  the  turbulent  case  and  thus 
are  given  in  reference  (1).  Similarly,  the  gas  phase 
mass  conservation,  ideal  gas  equation  as  well  as  the 
equations  for  determining  p  (assumed  spatially 
uniform)  and  u  remain  unchanged  and  are  not  repeated 
here.  Also,  the  procedure  to  determine  the  droplet 
surface  temperature  is  exactly  the  same  as  described 
in  reference  (1). 

To  complete  the  above  set  of  equations,  transport 
equations  for  k  and  e  are  required.  The  derivation 
of  these  equations  along  with  the  assumptions  in¬ 
volved  in  modeling  the  various  correlation  terms  will 
not  be  discussed  here  (see  references  6  and  10) .  The 
final  form  of  the  equations  for  one -dimensional  mean 


Several  investigators  have  studied  the  phenomenon  of 
particle  diffusion.  One  of  the  earliest  studies  is 
due  to  Tchen(12).  Under  highly  restrictive  condi¬ 
tions  which  are  not  applicable  for  a  real  case  (see 
references  13  and  14  for  details),  his  study  indi¬ 
cated  that  the  diffusion  coefficients  for  the  particle 
and  for  the  gas  are  the  same.  Further  studies  (IS, 16) 
showed  that  the  ratio  of  particle  diffusivity  to  eddy 
diffusivity  is  a  function  of  particle  response  time 
and  both  the  Lagrangian  and  Eularian  microscales  of 
turbulence.  Here  it  is  hypothesized  that  the  ratio 
of  the  gas  phase  diffusivity  and  the  particle  diffu¬ 
sivity  Dp  is  solely  dependent  on  the  ratio  of  a 
parti cle  response  time  and  a  turbulent  eddy  time, 
i .  e.  , 
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where  tpt  -  tp/tt. 

It  is  expected  that  if  tpt  0,  particle  diffu¬ 
sivity  will  be  equal  to  the  gas  phase  diffusivity. 

On  the  other  hand  for  a  very  large  tpt,  the  particle 
diffusivity  may  be  negligible.  The  particle  response 
time  may  be  expressed  in  terms  of  a  given  velocity 
perturbation  and  the  particle  acceleration,  i.e., 
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using  Stoke' s  drag  law  for  a  spherical  particle, 
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t  can  be  expressed  as. 
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Turbulence  is  generally  characterized  by  a  wide 
range  of  length  and  time  scales.  Large  scale  turbu¬ 
lence  may  be  characterized  by  the  length  and  time 
scales  associated  with  large  size  eddies  which  carry 
most  of  the  energy.  For  these  eddies,  the  character¬ 
istic  length  and  time  scales  may  be  defined  as  £t  = 
k5/2/e  and  tt  *  k/e.  Small  scale  turbulence  is 
represented  by  Kolmogorov's  microscales (12)  where 
the  time  scale  is  given  as  tt  *  (v/e)1/2.  in  order 
to  obtain  a  lower  bound  on  tpt,  it  is  appropriate  to 
consider  the  larger  of  the  two  turbulent  time  scales. 
Then  using  Eq.  [16],  t^t  can  be  written  as 
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which  can  be  rewritten  as 
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According  to  Eq.  (17),  the  particle  diffusivity 
depends  on  (a)  particle  mass  density/gas  phase  mass 
density,  (b)  square  of  the  ratio  of  particle  size 
and  a  turbulent  eddy  size  and  (c)  turbulent  diffusi¬ 
vity/  laminar  diffusivity.  Generally,  the  liquid 
density  is  two  to  three  orders  of  magnitude  higher 
than  the  gas  density.  If  the  turbulent  diffusivity 
is  one  to  two  orders  of  magnitude  higher  than  the 
laminar  one,  then  for  droplets  of  not  too  small  a 
size  (say  for  droplet  diameter  of  more  than  40 
microns),  tpt  will  usually  be  quite  large.  Under 
those  conditions,  one  can  assume  that  the  particle 
diffusion  is  negligible.  It  is  worth  noticing  that 
tpt  ■*  "^corresponds  to  the  case  considered  by 
Tchen ( * 2) .  It  may  also  be  mentioned  that  the  inclu¬ 
sion  of  particle  diffusivity  will  tend  to  give  a 
parabolic  character  to  the  liquid  phase  equations 
which  are  hyperbolic,  otherwise.  The  task  of  inclu¬ 
ding  the  particle  diffusion  is  being  pursued  by  the 
present  authors.  At  the  moment,  the  particle  diffu¬ 
sivity  is  neglected.  Then  the  liquid  phase  equations 
are  the  same  as  given  in  reference  (1). 

SOLUTION  PROCEDURE 

The  solution  procedure  used  here  is  more  or  less 
the  same  as  described  in  reference  (1).  The  non- 
dimensionalization  and  the  linearization  of  the 
governing  equations  are  obtained  in  the  same  manner. 
The  appropriate  boundary  conditions  for  the  gas  phase 
and  the  liquid  phase  variables  are  also  described  in 
reference  (1).  For  k  and  e  equations,  the  following 
boundary  conditions  are  used: 
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near  the  combustor  ends 


It  must  be  noted  that  the  above  relations  are  applied 
"near"  but  not  "at"  the  wall. 

The  linearized  equations  are  replaced  by  finite- 
difference  approximations  by  employing  an  implicit 
scheme,  similar  to  Crank-Nicolson  scheme.  This  leads 
to  a  tridiagonal  matrix  for  each  of  the  gas  phase 
variables  Yf,  Y0,  Yn,  K,  kf  and  e.  For  the  liquid 
phase  equations,  the  first  spatial  derivative  is 
approximated  as 
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where  i  and  j  represent,  respectively,  the  iAx  and 
j At  locations  in  space  and  time.  This  yields  an 
explicit  solution  of  the  liquid  phase  variables. 

The  initial  conditions  for  the  gas  phase  as  well  as 
for  the  liquid  phase  variables  are  provided  in  an  ad 
hoc  but  self-consistent  manner.  The  flame  initiation 
is  modeled  by  imposing  an  increase  in  gas  temperature 
and  a  decrease  in  fuel  and  oxidizer  concentrations  in 
a  small  zone  at  the  left  end  of  the  combustor.  The 
details  of  this  model  are  given  in  reference  (1). 

RESULTS 

Calculations  have  been  made  for  turbulent  flame 
propagation  in  homogeneous,  as  well  as  heterogeneous, 
mixtures  of  air  and  normal  decane  fuel.  In  both  the 
cases,  the  initial  conditions  are  identical  except 
that  in  the  homogeneous  case,  fuel  is  present  in  the 
vapor  form.  The  initial  pressure  is  one  atmosphere, 
and  the  overall  equivalence  ratio  is  0.95.  Several 
computer  runs  have  been  made  with  initial  k  and  c 
values  being  different  in  each  run.  The  variables 
are  plotted  versus  spatial  position  for  different 
times  to  indicate  the  propagation  of  flame.  The 
profiles  at  t  *  0  represent  the  initial  conditions 

Results  of  two  runs  (Pj  and  Pj)  for  the  turbu¬ 
lent  flame  in  homogeneous  mixtures  are  shown  in 
Figs.  1  and  2.  The  only  difference  in  the  two  runs 
is  in  the  initial  values  of  k  and  e  (see  Table  2). 
Profiles  of  temperature,  fuel,  and  oxidizer  mass 
fractions  in  Figs,  la  to  lc  clearly  indicate  a 
propagating  flame.  Figure  Id  gives  the  gas  velocity 
profile  at  t  *  60  ms .  The  gas  expansion  caused  by 
heat  release  as  a  result  of  chemical  reaction  is 
quite  evident  there.  One  of  the  more  notable  features 
of  runs  Pj  and  ?2  is  demonstrated  by  Figs,  le.  If, 
and  2b.  As  indicated  by  these  figures  and  further 
substantiated  by  Table  2,  the  profiles  of  k  and  c  at 
later  times  are  very  sensitive  to  the  initial  k/c 
ratio  which  represents  a  turbulence  decay  time.  This 
may  be  explained  by  examining  the  production  and 
dissipation  terms  in  k  (or  e)  equation.  If  the 
laminar  diffusivity  is  neglected,  using  Eq.  (10), 
the  source  terms  for  the  turbulent  energy  for  the 
pre-mixed  case  (m^  =  0)  can  be  written  as 

sk  .  Pk  -  e  (19) 

where  P  the  production  term,  is 

pk  -  Ikl7  (0-18lir-  11  (20) 

Since  3u/3x  term  is  mainly  determined  by  the  heat 
release  due  to  chemical  reaction,  it  is  significant 
only  in  the  flame  region.  As  a  result,  as  indicated 
by  Fig.  le,  k  increases  across  the  flame.  However, 
the  overall  k  value  throughout  the  combustor  shows  an 
initial  decrease.  This  is  due  to  the  fact  that  it 
takes  a  finite  time  (similar  to  ignition  delay  time) 
for  a  propagating  flame  to  become  establ ishea. 

During  that  time,  3u/3x  is  relatively  small  and  as  a 
result,  the  turbulence  decay  term  (e)  dominates  the 
turbulence  production  term.  After  that  initial 
period,  k  starts  increasing  with  time  (see  Fig.  le 
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and  Table  2).  The  turbulent  energy  dissipation  rate 
follows  the  sane  behavior. 
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TEMPERATURE  PROFILES 
Fig.  1o 


F*  1b 


Fig.  1c 


Fig  1* 


Fig  )f 


The  results  of  run  P2  tell  quite  a  different 
story.  As  indicated  by  Fig.  2b,  the  values  of  k  and 
e  not  only  decrease  monotonically  with  time,  but  they 
show  an  additional  dip  in  the  flame  region.  This 
occurs  since,  for  this  run,  the  production  term 
becomes  negative  in  the  flame  region.  According  to 
Eq.  (20),  the  sign  of  this  term  depends  on  the  rela¬ 
tive  amplitudes  of  k/e  (turbulence  decay  time)  and 
8u/3x,  related  to  the  inverse  of  the  chemical  reaction 
time.  For  example,  the  value  of  3u/3x  is  more  or 
less  the  same  for  runs  ?\  and  P2  as  the  chemical 
kinetics  is  identical  for  these  runs.  However,  the 
initial  k/e  in  run  P2  was  four  times  smaller  than 
that  in  run  Pi.  As  a  result,  k  increases  across  the 
flame  in  run  *1  and  decreases  across  the  flame  in  run 
P2.  The  dissipation  rate  (e)  follows  the  same 
behavior.  Hence,  the  profiles  of  k  and  e  are  mainly 
determined  by  the  turbulence  decay  time,  defined  by 
k/e,  and  the  turbulence  generation  time,  defined  by 
the  inverse  of  the  chemical  reaction  rate.  As  a 
result,  these  two  time  scales  greatly  influence  the 
turbulent  flame  propagation.  It  may  also  be  noted 
that  the  change  of  sign  for  the  production  term  in 
the  above  case  is  directly  attributed  to  the  inclu¬ 
sion  of  compressibility  effect  in  the  turbulent 
modeling . 


TEMPERATURE  PROFILES 
Fig.  2o 


Results  of  flame  propagation  in  a  heterogeneous 
mixture  are  given  in  Figs.  3  to  6.  Figure  3  gives 
the  results  for  laminar  as  well  as  turbulent  flame 
propagation.  It  may  be  mentioned  that  the  laminar 
results,  reported  here,  are  obtained  from  the  present 
turbulent  code  by  taking  turbulent  diffusivity  to  be 
zero.  These  results  are  identical  to  those  reported 
in  reference  (1).  It  is  noteworthy  that  inspite  of 
the  identical  initial  conditions,  the  flame  initia¬ 
tion  process  is  different  in  the  laminar  and  turbulent 
cases.  The  flame  initiation  time  (similar  to  an 
ignition  delay  time)  is  controlled  by  chemical  and 
diffusion  processes  (convection  is  absent  initially). 
Because  of  the  higher  diffusion,  it  takes  longer  for 
a  propagating  flame  to  become  established  in  the 
turbulent  case  (see  gas  temperature  profiles  at 
t  *  S  is  in  Fig.  3a).  Once  it  starts  propagating, 
the  turbulent  flame  moves  faster  than  the  laminar 
one.  The  expected  features,  like  increased  flame 
speed  and  thickness  for  the  turbulent  case,  are  quite 
evident  in  Figs.  3a  to  3e.  The  typical  values  are 


Diffusivity  in 
the  flame  region 
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Flame  Thickness 


LAMINAR 
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1.5  cm  /sec 
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80  cm/sec 
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NUMBER  DENSITY  profiles 
Fi»3« 


For  both  cases,  Fig.  3b  indicates  two  peaks  in 
the  fuel  vapor  mass  fraction  profiles.  The  first 
peak  represents  a  fuel  high  zone  left  behind  a  pre¬ 
mixed  type  propagating  flame,  represented  by  the 
second  case.  Because  of  the  higher  diffusion,  the 
peaks  for  the  turbulent  case  are  not  as  sharp  and  as 
high  as  for  the  laminar  case.  Also,  the  phenomenon 
of  a  secondary  diffusion  flame  (see  reference  1) 
seems  to  be  absent  for  the  turbulent  case.  For  the 
laminar  case,  the  secondary  flame  appears  after  t  = 
100  ms  and  is  not  shown  here.  The  profiles  of 
liquid  phase  properties  are  shown  in  Figs.  3d  and 
3e.  Two  points  are  worth  mentioning:  (i)  In  spite 
of  the  fact  that  the  particle  diffusivity  is  assumed 
to  be  zero,  the  liquid  phase  properties  show,  rela¬ 
tively,  a  smooth  variation  for  the  turbulent  case, 
(ii)  The  higher  droplet  surface  temperature  in  the 
hot  region  of  the  combustor  for  the  turbulent  case 
is  due  to  the  higher  pressure,  since  at  a  given  time, 
more  mass  has  been  burnt  for  the  turbulent  case. 

The  results  of  three  more  runs  for  the  turbulent 
case  are  shown  in  Figs.  4  to  6.  All  the  four  turbu¬ 
lent  runs  (runs  Tj,  T2,  T3,  and  T^)  differ  from  each 
other  only  in  their  initial  k  and  e  values.  Other 
initial  conditions  as  well  as  the  chemical  kinetics 
are  identical.  The  notable  feature  of  these  four 
runs  is  that  the  turbulent  flame  propagation  is 
greatly  influenced  by  the  initial  k/c  values. 
Depending  on  the  initial  k/e  value  (and  the  chemical 
kinetics) ,  one  can  have  an  accelerating  or  decelara- 
ting  flame.  Coupled  to  this  is  the  fact  that 
depending  on  k/e  and  chemical  kinetics  (through 
3u/3x  term),  the  turbulence  can  increase  or  decrease 
across  the  flame.  This  is  substantiated  by  the 
values  in  Table  2,  where  the  maximum  values  of  k  and 
e,  normalized  by  their  respective  initial  values,  as 
functions  of  time  are  tabulated.  Moreover,  as  these 
values  indicate,  the  time  rate  of  increase  (or 
decrease)  of  k  and  c  is  a  function  of  the  initial  k/e 
value.  For  run  T4,  initial  k/e  value  is  12.5  as 
compared  to  a  value  of  50  for  run  T2.  The  values  of 
k  and  e  decrease  relatively  faster  in  run  T4  as  com¬ 
pared  to  run  T2.  It  may  also  be  noted  that  for  run 
Tj,  the  maximum  values  of  k  and  e  first  decrease  and 
then  increase,  whereas  for  run  T3,  these  increase 
monotonically .  This  is  due  to  the  higher  initial 
diffusivity  for  run  Tj.  Because  of  higher  diffusion, 
it  takes  longer  for  a  flame  to  get  initiated  in  run 
Tj.  This  causes  the  turbulence  production  term  to  be 
initially  smaller  than  the  dissipation  term  for  run 
Tj.  A  similar  situation  has  been  mentioned  for  the 
pre-mixed  case. 
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FUEL  VAPOR  MASS  FRACTION  PROFILES 
Fig  5  b 


OXIDIZER  MASS  FRACTION  PROFILES 
Fig  5c 


Fig.  6« 

Some  additional  g-ross  features  of  the  turbulent 
flame  propagation  for  all  the  runs  are  summarized  in 
Table  3.  It  may  be  worth  noting  that  in  runs  T5  and 
T4,  the  flame  speed  and  thickness  are  roughly  the 
same.  This  is  in  spite  of  the  fact  that  initial 
diffusivity  in  run  T4  is  four  times  higher  than  that 
in  run  T3.  Again,  the  reason  is  that  the  turbulence 
level  increases  in  run  T,  and  decreases  in  run  T4. 

All  of  the  aforementioned  computer  runs  have 
been  made  on  DEC- 20  system.  In  order  to  limit  the 
computation  costs,  in  almost  all  the  runs,  computa¬ 
tion  has  been  terminated  when  flame  traverses  half 
the  combustor  length.  On  DEC-20  system,  this  takes 
typically  90  minutes  of  computer  time,  which  is 
equivalent  to  about  30  minutes  on  relatively  faster 
systems  (say  IBM  360).  In  all  the  runs,  the  step 
sizes  in  space  and  time  were,  respectively,  1 . 25  mm 
and  0.01  ms.  These  values  were  required  to  resolve 
the  combustor  process,  i.e.,  the  flame  thickness  and 
the  reaction  rate. 
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S UNWARY  AND  CONCLUSIONS 

The  main  findings  of  the  present  study  may  be 
summarized  as  follows: 

1.  The  laminar  flame  propagation  model  of  reference 
(1)  has  been  successfully  extended  to  predict 

the  gross  features  of  turbulent  flame  propagations 
in  air-fuel  drop-fuel  vapor  mixtures  in  a  one¬ 
dimensional  combustor.  The  turbulence  model  used 
is  based  on  the  gradient  transport  hypothesis, 
where  the  eddy  diffusivity  is  determined  by 
solving  the  transport  equations  for  k  and  c. 

2.  By  a  simple  consideration  of  a  particle  response 
time  and  a  turbulent  eddy  time,  an  estimate  is 
made  for  the  order  of  magnitude  of  particle  diffu¬ 
sivity.  It  is  shown  that  the  ratio  of  particle 
diffusivity  to  gas  phase  turbulent  diffusivity 
depends  on  (a)  particle  mass  density/gas  phase 
mass  density,  (b)  square  of  the  ratio  of  a  particle 
size  to  a  turbulent  eddy  size,  and  (c)  turbulent 
diffusivity/laminar  diffusivity.  Conditions  are 
discussed  where  the  particle  diffusivity  may  be 
neglected.  It  is  neglected  presently  but  can  be 
included  in  a  future  study. 

3.  The  characteristics  of  turbulent  flame  propagation 
in  both  homogeneous  and  heterogeneous  mixtures 
seem  to  be  highly  sensitive  to  the  relative 
magnitudes  of  k/e,  which  defines  a  turbulence 
decay  time,  and  inverse  of  the  chemical  reaction 
rate,  which  defines  a  turbulence  generation  time. 
This  conclusion  should  be  examined  more  generally 
by  further  numerical  studies  and  experimental 
results . 

4.  The  present  model  has  been  used  for  a  mono-dis- 
persed  spray.  The  poly-disperse  spray  can  be 
considered  in  a  future  study. 
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Effect  of  Numerical  Modeling  on 
One-Dimensional  Enclosed 
Homogeneous  and  Heterogeneous 
Deflagrations 

The  objective  of  this  study  is  to  evaluate  the  effectiveness  of  several  finite  difference 
techniques  in  modeling  unsteady  homogeneous  and  heterogeneous  flames  in  closed 
combustors.  The  problem  of  propagation  of  a  pre-mixed  homogeneous  flame  and  a 
heterogeneous  flame  in  air-fuel  droplet-fuel  vapor  mixture  is  formulated.  The 
governing  equations  have  been  solved  by  using  four  different  numerical  techniques. 
For  the  pre-mixed  case,  flame  propagation  is  computed  for  two  fuels,  propane  and 
n-octane.  For  the  two-phase  heterogeneous  case,  n-octane  fuel  is  used.  The  results 
are  obtained  in  terms  of  the  profiles  of  gas  and  liquid  phase  properties  at  different 
instances  of  time  as  the  flame  traverses  the  combustor.  From  these  results,  the 
relative  efficiency  and  accuracy  of  each  of  the  numerical  techniques  is  discussed. 


INTRODUCTION 

Numerical  computation  of  unsteady  flows  with 
combustion,  where  the  flow  properties  are  assumed  to 
vary  in  one  or  two  spatial  directions,  have  reached 
an  advanced  level  of  sophistication.  Both  explicit 
and  implicit  methods  have  been  used  to  a  variety  of 
flow  situations.  Generally,  the  implicit  numerical 
methods  have  been  reconmended  because  of  their  favor¬ 
able  stability  properties.  The  explicit  schemes  for 
a  parabolic  unsteady  PDE  are  required  to  satisfy  cer¬ 
tain  stability  requirements,  which  impose  a  restric¬ 
tion  on  the  temporal  step  size  relative  to  the  spatial 
mesh  size.  In  one-dimensional  problems,  these  sta¬ 
bility  conditions  are  expressed  at  At<^x/|u|+c  and 
t  i  ax2/ 2i).  Thus,  the  maximum  allowable  time  step 
could  be  dictated  by  the  spatial  mesh  size,  rather 
than  by  a  characteristic  time  of  a  physical  process 
under  consideration.  Since  the  implicit  methods  are 
generally*  free  from  these  restrictions,  it  may  be 
possible  to  use  a  relatively  larger  time  step  than 
what  Is  permissible  in  explicit  schemes.  However, 
the  advantage  of  using  an  implicit  formulation  may 
prove  illusory  in  many  practical  situations.  An 
implicit  scheme  requires  a  relatively  larger  number 
of  arithmetic  operations  than  an  explicit  scheme. 

For  example,  as  discussed  by  Cheng  [3],  the  solution 
of  a  scalar  PDE  will  require  five  times  more  compu¬ 
tational  effort  per  time  step  in  an  implicit  formu¬ 
lation  than  In  an  explicit  one.  For  a  coupled  non- 

*StfTlY  Implicit  schemes  are  required  to  satisfy 
a  cell  Reynolds  number  criterion  and  any  stability 
conditions  arising  from  the  implementation  of  boundary 
conditions. 
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linear  PDE  set,  the  implicit  formulation  will  be  even 
more  time  consuming  because  an  iterative  scheme  is 
required  to  take  care  of  the  non-linearities  and  the 
coupling  in  the  governing  equations.  As  a  result,  an 
implicit  method  can  retain  its  computational  advantage 
only  if  the  time  step  is  considerably  larger  (say  an 
order  of  magnitude  and  higher)  than  that  allowed  by 
stability  requirements. 

In  an  effort  to  make  the  implicit  formulation 
more  efficient,  several  linearization  techniques  have 
been  developed.  These  techniques  mainly  fall  into 
two  categories.  One  is  the  quasi-linearization 
(Newton-Raphson)  iterative  process  as  originally  sug¬ 
gested  by  Bellmann  and  Kalaba  [1],  and  later  used  by 
Sharma  and  Sirignano  [4]  for  solving  a  non-linear  POE 
set.  In  this  method,  which  requires  the  solution  of 
a  tridiagonal  matrix  for  each  of  the  scalar  equations, 
the  linearization  is  conducted  by  using  a  Taylor 
expansion  in  an  iterative  space.  The  second  class  of 
linearization  processes  has  been  suggested  indepen¬ 
dently  by  Briley  and  McDonald  [2]  and  Lindenmuth  and 
Killer  [5].  The  technique,  which  has  been  called  a 
time-1 inearization,  uses  a  Taylor  expansion  in  time 
and  requires  the  solution  of  a  block  tri-diagonal 
system  as  opposed  to  a  simple  tri-diagonal  system  in 
quasi-linearization  scheme.  The  main  attraction  of 
the  time- 1 ineari zation  method  over  the  quasi -1 ineari - 
zation  method  is  that  the  former  leads  to  a  non-iter¬ 
ative  algorithm.  However,  for  any  of  these  implicit 
1 inearization  methods  to  have  an  edge  over  the  expli¬ 
cit  schemes,  the  temporal  step-size  still  has  to  be 
at  least  an  order  of  magnitude  larger  than  the  maximum 
permissible  in  the  explicit  schemes.  On  the  other 
hand,  there  are  many  practical  situations  where  the 
space  and  time  step  restrictions  imposed  due  to  the 
finite  differencing  of  the  governing  equations  are  not 
as  important  as  those  imposed  by  the  physical  and 
chemical  processes  under  consideration.  In  such  cases, 
the  advantages  of  even  highly  efficient  implicit 
schemes  become  questionable.  The  physical  situation 
considered  in  the  present  study  falls  into  this 


category. 

An  unsteady  propagation  of  pre-mi xed  homogeneous 
flames  as  well  as  of  heterogeneous  flames  in  an  air- 
fuel  droplet-fuel  vapor  mixture  in  enclosed,  one¬ 
dimensional  combustors  is  formulated.  For  the  pre¬ 
mixed  homogeneous  case,  the  set  of  parabolic  PDE‘s 
is  solved  by  four  different  finite-difference  schemes; 
an  explicit  approach,  an  iterative  sequential  implicit 
method  without  any  formal  linearization,  an  iterative 
quasi-linearization  implicit  method  with  a  scalar  tri¬ 
diagonal  matrix  system  and  a  non-iterative  time- 
linearized  implicit  technique  with  a  block  tri-diago¬ 
nal  matrix.  In  each  of  the  four  methods,  a  central 
or  upwind  difference  scheme  is  used  for  the  spatial 
derivatives,  depending  upon  the  local  cell  Reynolds 
number.  The  relative  merits  of  each  of  the  numerical 
techniques  is  discussed. 

For  the  case  of  heterogeneous  flame  propagation 
in  a  two-phase  flow,  the  gas-phase  equations,  a  sub¬ 
set  of  parabolic  PDE*s,  is  again  solved  by  the  same 
four  numerical  techniques.  The  liquid  phase  equa¬ 
tions,  a  subset  of  hyperbolic  PDE’s,  is  solved 
explicitly.  It  is,  however,  worth  mentioning  that 
because  of  the  coupling  between  the  two  sets,  all  the 
implicit  schemes  require  iteration  and,  in  each  iter¬ 
ation,  the  liquid  phase  properties,  although  obtained 
explicitly,  are  updated.  In  the  next  section,  the 
physical  model  and  governing  equations  are  presented. 
This  is  followed  by  the  discussion  of  the  solution 
procedure  employed  in  each  of  the  methods.  Finally, 
the  results  are  presented  and  the  relative  merits  of 
the  various  numerical  techniques  are  discussed. 

PHYSICAL  MODEL  AND  GOVERNING  EQUATIONS 


The  physical  situation  considered  in  the  present 
study  is  that  of  a  one-dimensional  combustor,  enclosed 
at  both  ends  for  both  the  pre-mi xed  and  the  hetero¬ 
geneous  two-phase  cases.  The  initial  profiles  of  the 
gas  phase  properties,  as  well  as  of  the  liquid  phase 
properties,  are  prescribed  (see  Figures  1-10).  The 
conservation  of  gas  phase  properties  is  described  by 
the  following  equations: 


ao/at  +  a/ ax  (pu)  *  Sp  (1 ) 

p  aY./at  ♦  pu  aY^/ax  -  a/3x  (u  aY^/ax)  *  si  (2) 
os$/at  +  pu  a$/ax  -  a/ax  (ya^/ax)  =  (3) 
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where  p,  u,  T,  p  and  Y.  are  respectively  the  density, 
velocity,  temperature,  pressure  and  mass  fractions 
for  the  gaseous  phase.  Furthermore,  y,  R,  y  and 
are  the  gas  viscosity,  universal  gas  constant,  ratio 
of  specific  heats,  and  molecular  weight  for  species  i. 
Also,  x  and  t  are  the  spatial  and  temporal  coordinates 
and  L  is  the  combustor  length.  Sp,  S-j ,  and  S$  are  the 
source  or  sink  terms  for  the  gas  phase  and  are  pre¬ 
scribed  in  Table  1.  These  are  discussed  in  greater 
detail  in  [6].  In  deriving  the  above  equations,  it 
is  assumed  that: 

1.  the  volume  occupied  by  the  liquid  phase  is 
negl igible, 

2.  pressure  is  only  a  function  of  time, 

3.  the  species  diffusion  follows  Fick's  law 
with  equal  mass  di ffusi vi ties  for  all  species, 

4.  Prandtl  and  Schmidt  numbers  are  unity, 

5.  radiative  heat  transfer  is  negligible, 

6.  viscous  dissipation  is  negligible  and  the 
gas  kinetic  energy  is  small  as  compared  to 
its  thermal  energy, 

7.  gas  phase  is  thermally  and  calorically  per¬ 
fect  with  equal,  constant  specific  heat,  and 

8.  the  gas  phase  chemical  reaction  is  governed 
by  a  single-step  second-order  kinetics. 

LIQUID-PHASE  EQUATIONS 


The  major  assumption  in  describing  the  liquid 
phase  is  that  the  droplet  properties  are  described 
by  an  Eulerian  approach.  Further,  a  single  component 
mono-disperse  dilute  fuel  spray  is  considered.  As  a 
result,  the  interactions  between  the  droplets  are 
neglected.  The  droplets  are  assumed  to  be  spherical 
with  the  convection  effect  being  described  by  a  Ranz- 
Marshall  correlation  [7].  Furthermore,  the  drop^t 
interior  temperature  is  considered  uniform.  This  is 
equivalent  to  a  large  liquid  thermal  conductivity 
assumption.  The  derivation  of  source  (or  sink)  terms 
for  the  gas  and  liquid  phases  are  based  on  this 
assumption.  Then,  the  subset  of  hyperbolic  PDE's  for 
the  liquid-phase  properties  are: 
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where  n,  ut,  Re,  and  Te  are  respectively,  the  droplet 
number  density,  velocity,  radius,  and  surface  temper¬ 
ature.  Furthermore,  o*  and  Cp.  are  the  density  and 
specific  heat  for  the  liquid  fuel.  The  equilibrium 
vapor  fuel  mass  fraction  { Y ^  )  as  a  function  of  T-  is 
prescribed  by  Antoine's  relation  (13,14)  as  discussed 
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in  [8].  Finally,  the  drag  law  (12)  is  used  in  a  form 
suggested  by  Ingebo  [9]. 


where  *  stands  for  any  of  the  variables  and  subscript 
x  denotes  the  spatial  derivative. 


It  is  noteworthy  that  in  spray  combustion 
studies,  the  controversy  regarding  the  existence  of 
envelope  flame  surrounding  the  individual  droplets 
still  exists.  The  experimental  studies  by  Chigier 
[10]  show  no  evidence  of  individual  droplet  flames. 
However,  under  certain  conditions  such  as  large  drop¬ 
lets  in  hot  surroundings,  the  existence  of  an  envelope 
flame  is  a  distinct  possibility.  In  the  present  study 
both  of  the  cases  are  considered,  where  the  presence 
of  a  droplet  envelope  flame  is  determined  by  a  cri¬ 
tical  Damkohler  number  as  used  by  Law  [11]. 

NUMERICAL  CONSIDERATIONS  AND  SOLUTION  PROCEDURES 

Here  we  describe  the  method  of  solution  for  each 
of  the  numerical  schemes.  The  computer  codes  have 
been  developed  in  such  a  way  that  the  same  code  can 
be  used  for  the  pre-mixed  single-phase  and  hetero¬ 
geneous  two- phase  flames.  The  pre-mixed  case  is  re¬ 
covered  from  the  two  phase  case  by  changing  the  ini¬ 
tial  conditions  and  making  S0  identically  equal  to 
zero.  Therefore,  the  solution  procedure  is  outlined 
only  for  the  latter  case.  In  all  the  numerical 
schemes,  the  general  procedure  is  as  follows: 

First  of  all,  the  liquid-phase  properties  are 
advanced  in  time  explicitly  by  using  the  gas  and 
liquid  phase  properties  at  the  known  time  level. 

Since  the  liquid-phase  equations  are  hyperbolic  in 
nature,  the  spatial  derivatives  are  finite  differenced 
by  using  an  upwind  scheme.  The  advanced-liquid  phase 
properties  are  then  used  to  integrate  the  gas-phase 
energy  and  species  conservation  equations.  In  all 
the  numerical  schemes,  the  convection  term  in  the  gas- 
phase  equations  is  f ini te-differenced  by  using  a  cen¬ 
tral  or  upwind  scheme  as  determined  by  the  local  cell 
Reynolds  number.  At  locations  where  the  cell  Reynolds 
number  exceeds  two,  the  difference  approximation  is 
changed  from  central  to  upwind  [12].  The  values  of 
$  and  Y .  obtained  at  the  advanced  time  level  are  then 
employed  in  (5,6)  to  obtain  new  gas  pressure  and  den¬ 
sity.  Finally,  Eq.  (7)  is  used  to  give  the  new  gas 
velocity  and  the  whole  cycle  is  repeated  to  advance 
to  the  next  time-step. 

The  gas-phase  equations  are  solved  by  using  four 
different  numerical  schemes.  In  the  explicit  schemes, 
the  solution  is  advanced  by  using  the  spatial  differ¬ 
encing  at  the  explicit  time- level.  In  the  fully 
implicit  scheme,  there  is  no  formal  linearization. 

The  finite  differencing  of  the  spatial  derivative  is 
written  at  the  implicit  time  level.  The  finite- 
difference  form  is  assumed  to  be  linear  and  the  solu¬ 
tion  is  obtained  by  inverting  the  tri-diagonal  matrix. 
The  equations  are  solved  in  a  sequential  manner. 

First  of  all,  the  energy  equation  is  solved,  which 
provides  an  updated  value  of  $,  which  in  turn  is  used 
for  solving  the  equation  for  Y..  Then  updated  values 
of  $  and  Y,  are  used  to  solve  for  YQ.  This  procedure 
of  sequential  updating  is  continued uin  an  iterative 
manner  until  the  desired  convergence  is  obtained  and 
the  cycle  is  repeated  to  advance  the  solution  to  the 
next  time-step. 

In  the  implicit  method  with  quasi-linearization 
in  an  iterative  space,  the  linearization  procedure 
can  be  described  by  considering  any  of  the  governing 
equations  as: 

n/n  *  zU,  *x,  *xx)  (15) 


Equation  (15)  is  linearized  as: 

(3*/3t)n+1  =  2°  +  (3Z/3*)n  Ati,n+1 
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Finite  differencing  of  (16),  along  with  the  boun¬ 
dary  conditions  (no  heat  and  mass  diffusion  at  the 
combustor  boundaries)  yields  a  tridiagonal  matrix  for 
each  of  the  variables.  The  method  still  requires 
Iteration  because  the  governing  equations  are  coupled. 
The  idea  of  using  the  quasi-linearization  process  over 
the  straight  forward  implicit  scheme  without  any  for¬ 
mal  linearization  is  to  expedite  the  convergence  pro¬ 
cess. 


The  time-linearization  non-iterative  formulation 
is  described  in  detail  by  Briley  and  McDonald  [2]. 

Only  a  brief  description  of  the  technique  is  given 
here.  For  a  first  order  non-linear  equation 

(3Y/3t)n+1  =  F(Y)  3/3x  G(Y)  (18) 

the  Taylor  series  expansion  of  non-linear  terms  in 
time  is  written  as 

(3Y/3t)n+1  =  (FGx)n 

*  [Fn(3Gx/3Y)n  +  Gx"  (3F/3Y)n] 

(Y^1  -  Vin)/At  (19) 

The  non-linear  source  terms  and  the  second  order  spa¬ 
tial  derivative  terms  are  treated  in  the  same  way. 
Using  an  appropriate  finite  differencing  and  the  boun¬ 
dary  condition  leads  to  a  tri diagonal  matrix.  For  a 
coupled  set  of  non-linear  equations  (i.e.,  when  Y  is 
a  column  vector),  the  above  linearization  process  will 
result  in  a  block  tri-diagonal  system,  which  can  be 
solved  either  by  a  standard  block  tri-diagonal  system 
solver  or  by  a  band  matrix  solver.  The  block  tri- 
diagonal  algorithm  is  similar  to  that  used  for  solving 
a  scalar  tri-diagonal  system  except  that  each  element 
itself  becomes  a  square  matrix.  The  band  matrix  sol¬ 
ver  is  a  more  efficient  way  of  inverting  a  block  tri¬ 
diagonal  matrix,  where  all  the  diagonals  in  the  block 
tri-diagonal  system  are  transformed  to  form  another 
matrix  and  the  solution  is  obtained  by  inverting  the 
transformed  matrix.  It  may  be  mentioned  that  in  the 
present  calculations,  the  band  matrix  solver  was  fifty 
percent  more  efficient  than  the  block  tri-diagonal 
matrix  solver. 

RESULTS  AND  COMPARISONS  OF  NUMERICAL  SCHEMES 

The  results  for  the  heterogeneous  flame  propaga¬ 
tions  in  air-fuel  vapor-fuel  droplet  (n-octane)  mix¬ 
ture  in  a  one-dimensional  enclosure  combustor  of 
length  10cm  are  presented  in  Figures  1-10.  To  keep 
the  computer  cost s  under  control,  computations  have 
been  made  until  the  flame  traverses  about  half  the 
combustor  length.  In  addition,  for  the  purpose  of 
clarity,  the  profiles  of  the  various  gas  and  liquid- 
phase  properties  are  plotted  for  only  half  the  length. 
The  initial  profiles  (at  t=0)  of  all  the  gas  and 
liquid-phase  properties  are  assumed  uniform  except 
that  a  flame  initiation  or  ignition  source  of  length 
0.5cm  and  of  temperature  1500°K  is  provided  near  the 
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left  end  of  the  combustor.  The  initial  value  of  all 
variables  are  varied  linearly  from  the  hot  zone 
values  to  ambient  values  in  the  next  0.5cm  to  avoid 
large  gradients  near  the  ignition  source.  The  ini¬ 
tial  pressure  is  one  atmosphere  and  the  overall  equi¬ 
valence  ratio  is  .98.  The  spatial  and  temporal  sizes 
(ix,At)  for  most  of  the  calculations  are  respectively 
,125cm  and  10*5  seconds.  The  spatial  step  size  (ax) 
was  selected  so  that  the  results  are  relatively  insen¬ 
sitive  to  a  further  reduction  in  ax.  The  largest 
permissible  At  for  the  present  computation  is  con¬ 
trolled  by  the  reaction  time  scale,  which  is  about 
0.1ms  for  n-octane  fuel.  It  is  important  to  point 
out  that  for  these  step  sizes,  all  four  aforementioned 
numerical  schemes  produce  results,  which  are  indis¬ 
tinguishable  from  each  other;  the  maximum  gas  temper¬ 
ature  difference  between  the  numerical  methods  is  less 
than  two  percent. 

The  gas  temperature  profiles  are  shown  in  Figure 
1.  As  this  Figure  indicates,  the  flame  initiation 
or  ignition  process  is  rather  slow,  taking  more  than 
20ms.  However,  after  this  initiation  or  ignition 
period,  the  flame  propagates  at  a  much  faster  rate. 

It  is  noteworthy  that  in  the  hot  region  behind  the 
propagating  flame,  the  lowest  gas  temperature  is  near 
the  wall.  This  is  due  to  the  nature  of  the  initial 
temperature  profile  used  for  ignition.  From  these 
temperature  profiles,  it  seems  rather  difficult  to 
distinguish  between  a  pre-mi xed  homogeneous  and  two- 
phase  heterogeneous  cases.  This  distinction  is  shown 
more  clearly  by  the  fuel  vapor  mass  fraction  profiles, 
plotted  in  Figure  2.  Because  of  the  slow  flame  ini¬ 
tiation  process  and  rather  fast  vaporization  process 
for  n-octane  fuel  (with  relatively  high  volatility), 
a  substantial  amount  of  fuel  vapor  is  left  behind  the 
propagating  flame.  This  is  clearly  exhibited  by  the 
first  peak  in  the  fuel  vapor  mass  fraction  profiles. 
For  a  relatively  less  volatile  fuel,  the  vaporization 
process  will  be  slow  and  we  will  not  get  such  a  high 
peak.  The  same  effect  could  also  be  obtained  by 
making  the  flame  initiation  process  faster.  It  should 
also  be  noted  that  the  peak  value  of  Y^  in  these 
profiles  continues  to  increase  upto  t  *  60ms.  By  that 
time,  the  droplets  in  the  hot  region  are  completely 
vaporized  (see  Figure  5)  and  consequently  the  peak 
value  of  Yj-  starts  decreasing  due  to  diffusion.  The 
second  peak,  which  is  more  clearly  visible  at  t  = 

60ms  and  80ms,  occurs  at  the  propagating  flame.  The 
oxidizer  mass  fraction  profiles  in  Figure  3  portray 
the  propagating  nature  of  the  flame  more  distinctly. 
Since  there  is  no  oxidizer  left  behind  the  flame,  the 
amount  of  fuel  vapor  in  that  region  remains  more  or 
less  unchanged.  The  gas  velocity  profiles  are  given 
in  Figure  4.  Across  the  propagating  flame  region,  the 
gas  velocity  has  the  largest  value,  (HOcm/sec)  and 
changes  sign.  This  is  also  the  region  where  the  gas 
and  liquid  phase  properties  have  the  steepest  gradi¬ 
ents,  and  it  is  due  to  the  resolution  of  this  region 
that  the  explicit  scheme  is  superior  to  the  implicit 
schemes.  Corresponding  to  this  maximum  velocity,  the 
CFL  (Courant-Frfedrichs-Lewy)  stability  condition 
requires  a  time  step  of  1ms,  which  is  the  most  strin¬ 
gent  stability  requirement;  the  diffusion  stability 
limit  being  At  larger  than  5ms.  Since  the  character¬ 
istic  time  of  the  chemical  process  in  the  present 
study  is  about  0.1ms,  the  time  step  is  completely 
dictated  by  this  process.  This  makes  the  explicit 
numerical  method  far  superior  than  any  of  the  implicit 
methods  used  here.  For  example,  for  the  same  amount 
of  computation,  the  computer  times  for  the  explicit 
scheme,  iterative  implicit  without  formal  lineariza¬ 


tion,  iterative  implicit  with  quasi-linearization,  and  9 

block  implicit  with  time-linearization  respecti vely , 
are  8,  32,  33,  and  45  minutes.  Thus,  the  explicit 
method  is  the  most  efficient  one  and  the  block  impli¬ 
cit  is  the  least  efficient  one.  The  inversion  of 
block  tridiagonal  matrix,  even  by  an  efficient  band 
matrix  solver,  is  a  highly  time-consuming  process  and 
it  makes  the  block  implicit  scheme  the  least  favor-  9 

able. 

Figure  5  gives  the  variation  of  droplet  radius 
with  spatial  locations  and  time.  The  valleys  in  these 
profiles  indicate  that  the  gas  temperature  is  the 
highest  in  these  regions.  After  t  -  60ms,  all  the 
droplets  in  the  hot  region  are  completely  vaporized. 

The  droplet  radius  profiles  at  t  =  60  and  80ms  indi-  • 

cate  an  interesting  phenomenon  at  the  spatial  location, 
where  the  profile  curvatures  change  sign.  For  the 
calculations  presented  in  Figs.  1-6,  a  droplet  igni¬ 
tion  criterion  [11]  has  been  used  to  determine  whether 
there  is  a  flame  surrounding  individual  droplets.  To 
the  left  of  this  spatial  location,  this  criterion  is 
satisfied  and  therefore  individual  droplets  have 
envelope  flames  in  that  region.  As  a  result,  the  ® 

vaporization  rate  is  also  relatively  high.  As  we  go 
to  the  right  of  this  spatial  location,  the  droplet 
ignition  criterion  is  no  longer  satisfied,  there  are 
no  envelope  flames,  and  consequently  the  vaporization 
rate  and  rate  of  decrease  of  droplet  size  reduces. 

Thus  the  proper  identification  of  the  droplet  ignition 
phenomenon  may  be  very  important  in  spray  combustion  ^ 

studies.  As  shown  in  Fiqure  6,  the  droplet  density 
number  increases  to  the  right  of  the  propagating 
flame  because  of  the  droplet  velocities,  induced  by 
the  gas  velocity.  However,  the  total  number  of  drop¬ 
lets  in  the  combustor  does  not  change.  The  droplet 
velocity  profiles  are  similar  to  the  gas  velocity 
profiles  and  are  not  shown  here. 

» 

Figure  7  shows  the  effect  of  changing  cold  gas 
temperature  on  the  two-phase  flame  propagation.  A 
comparison  of  Figures  1  and  7  indicates  that  increas¬ 
ing  the  cold  gas  temperature  from  300°K  to  400°K 
causes  a  marked  increase,  not  only  in  the  flame  ini¬ 
tiation  process,  but  also  in  the  flame  propagation 
speed.  Furthermore,  the  gas  temperatures  in  the  hot 
region  are  higher  in  Figure  7  than  in  Figure  1,  be-  • 

cause  the  gas  pressure  increases  at  a  faster  rate  in 
the  former  case. 

Figure  8  gives  the  temperature  for  the  case  when 
there  is  no  envelope  flame  surrounding  any  individual 
droplet.  This  case  is  obtained  by  simply  assuming 
that  the  droplet  ignition  criterion  is  never  satisfied  _ 

and  therefore  the  droplets  are  only  vaporizing.  A  9 

comparison  of  these  profiles  with  those  when  there 
are  individual  burning  droplets  (Figure  7)  indicates 
that  flame  propagation  is  faster  in  the  latter  case. 

A  simple  explanation  is  that  in  the  regions  where 
droplet  ignition  criterion  is  satisfied,  the  total 
burning  (combusting)  rate  is  much  larger  than  in  the 
non-burning  droplet  case  and  this  causes  an  increase  t 

in  the  flame  speed. 

Another  interesting  comparison  of  the  burning  and 
non-burning  droplet  cases  is  demonstrated  in  Figure  9, 
where  the  droplet  radii  are  plotted  versus  spatial 
locations  at  different  instances  of  time.  For  the 
completely  non-burning  droplet  case,  the  rate  of  de¬ 
crease  of  droplet  is  smaller  than  the  corresponding  9 

case  when  there  are  individual  burning  droplets.  This 
is  because  the  vaporization  rate  for  the  burning 
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droplets  is  much  higher  (about  twice  as  high  in  the 
present  study)  than  that  for  the  purely  vaporizing 
droplets.  Also,  the  change  in  curvature  of  the  drop¬ 
let  radius  profiles  as  discussed  earlier  for  the 
burning  droplet  case,  is  not  observed  when  the  drop¬ 
lets  are  purely  vaporizing.  The  conclusions  regarding 
the  efficiencies  of  numerical  schemes  remain  unaltered 
for  the  cases:  a)  when  cold  gas  temperature  is  raised 
to  400° K,  and  b)  when  droplets  are  considered  to  be 
only  vaporizing. 

The  aforementioned  numerical  schemes  have  also 
been  used  to  compute  the  homogeneous  pre- mixed  flame 
propagation.  Again  the  time-step  restrictions  are 
imposed  by  the  chemical  reaction  rather  than  by  any 
stability  limits.  To  examine  this  point  further, 
computations  are  made  for  n-octane  and  propane  fuels, 
which  have  different  chemical  kinetics  parameters. 

The  temperature  profiles  for  these  fuels  are  given 
in  Figure  10.  For  n-octane  fuel,  the  same  amount  of 
computation  for  the  explicit,  iterative  implicit 
without  formal  linearization,  iterative  implicit  with 
quasi-1 inearization,  and  block  implicit  schemes  re¬ 
quires  respectively  5,  25,  25.5  and  36  minutes  of 
computer  time.  For  propane  fuel,  the  superiority  of 
the  explicit  methods  increases.  This  is  because  the 
characteristic  reaction  time  is  smaller  for  propane 
than  for  n-octane.  As  a  consequence,  the  flame 
initiation  as  well  as  flame  propagation  are  faster 
for  propane  fuel.  In  addition,  the  gas  properties 
have  relatively  steeper  gradients  in  the  flame  region 
and  all  the  implicit  schemes  take  a  large  number  of 
iterations  to  compute  this  region  with  an  acceptable 
degree  of  accuracy. 

It  seems  appropriate  to  compare  some  of  the 
present  results  with  those  presented  in  an  earlier 
study  [6],  which  considered  the  flame  propagation  in 
air-  n-decane  fuel  droplet  mixture  with  non-uniform 
droplet  interior  temperature.  The  results  of  that 
study  indicated  the  existence  of  a  secondary  diffusion¬ 
like  flame  behind  the  primary  propagating  flame.  The 
secondary  flame  is  not  observed  in  the  present  calcu¬ 
lations.  The  reason  is  the  highly  volatile  nature 
of  n-octane,  which  always  creates  a  fuel  rich  zone  in 
the  vicinity  of  the  propagating  flame.  For  n-decane 
fuel,  which  is  not  so  volatile,  a  fuel  lean  zone  is 
formed  near  the  propagating  flame.  This  causes  some 
unbumt  oxidizer  to  be  left  behind  the  primary  flame, 
which  subsequently  diffuses  towards  the  initially 
formed  fuel  rich  zone  (near  the  left  end  of  the  com¬ 
bustor)  and  forms  a  secondary  diffusion  flame. 

CONCLUSIONS 

(1)  An  unsteady  flame  propagation  in  a  two-phase 
heterogeneous  mixture  is  studied  by  using 
four  different  finite-difference  schemes; 

an  explicit  method,  an  iterative  implicit 
method  without  any  formal  linearization,  an 
iterative  implicit  method  with  quasi-1 inear¬ 
ization,  and  a  time-1 inearization  block 
implicit  method. 

(2)  The  results  indicate  that  the  explicit 
scheme  is  by  far  the  best.  This  is  mainly 
due  to  the  fact  that  the  time  step  size  is 
controlled  by  the  chemical  reaction  and 
vaporization  rates,  and  not  by  any  stability 
requirements. 

(3)  Amongst  the  implicit  methods,  the  straight¬ 
forward  implicit  scheme  without  any  formal 


linearization  is  the  most  efficient  for  the 
present  calculations.  Thus,  a  quasi-linear¬ 
ization  process  does  not  expedite  the  conver¬ 
gence  of  the  solutions,  at  least  for  the 
problem  considered  in  the  present  study. 

The  block  implicit  scheme  is  the  most  expen¬ 
sive  one,  as  the  inversion  of  a  block  tri¬ 
diagonal  matrix  even  by  an  efficient  band- 
matrix  solver  is  a  time-consuming  process. 

(4)  Conclusions  (2)  and  (3)  should  not  be  ac¬ 
cepted  in  a  general  sense.  Implicit  schemes 
can  become  attractive  if  (a)  a  very  high 
spatial  resolution  is  needed  and/or  a  locally 
refined  mesh  is  used;  the  latter  situation 
arises  in  flows  with  more  than  one  length 
scale,  as  is  the  case  for  flows  with  thin 
boundary  layers,  (b)  the  temporal  accuracy 

is  of  no  importance;  for  example,  in  the 
computation  of  steady  state  flows  via  asymp¬ 
totic  temporal  approach,  and  (c)  the  fuel 
used  is  not  so  volatile.  The  block  implicit 
scheme  may  have  an  advantage  in  situations 
where  a  large  number  of  iterations  in  an 
iterative  implicit  scheme,  is  required. 
Furthermore,  the  block  implicit  scheme  should 
be  used  when  the  number  of  PDF's  in  the  block 
is  equal  to  the  number  of  flow  variables  to 
be  solved.  Such  was  not  the  case  in  the 
present  formulation  and  consequently,  an 
iteration  was  required  even  for  the  block 
scheme. 

(5)  The  physical  process  modeled  here  is  inher¬ 
ently  stiff,  e.g.,  the  various  time  scales 
are  quite  disparate.  In  this  situation,  the 
use  of  an  operator  splitting  should  be  exam¬ 
ined.  Also,  the  method  of  lines  should  be 
explored. 

(6)  Recently,  an  Eulerian-Lagrangian  scheme  has 
been  successfully  used  to  solve  a  two-phase, 
two-dimensional  jet  flow  [13].  It  is  recom¬ 
mended  to  explore  the  employment  of  such  a 
scheme  for  the  physical  situation  presented 
here. 
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TABLE  1  -  SOURCE  TERM  FOR  THE  GAS-PHASE  EQUATIONS 
(1-7) 
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Fig.l  Gas  Temperature  Profiles  for  Two-Phase  Case 


*0  =  w0  M0  *  S0V  where 

M0  =  5  for  propane 
e  12.5  for  octane 
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*  0  for  vaporizing  droplets 


Fig. 2  Fuel  Vapor  Mass  Fraction  Profiles,  Two-Phase 
Case 
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abstract 

A  multi-orifice  impulsed  spray  generator  with  two 
different  modes  of  operation  is  developed.  The  impulse 

iet  operation  mode  covers  a  wide  range  of  droplet 

frequencies,  and  produces  arrays  of  droplets  with 
controllable  spacing,  but  with  a  wide  sire  spectrum.  This 
ejector  also  operates  in  pressurized  flow  mode.  In  this 
mode,  optimum  disturbance  frequency  breaks  up  the  liquid 
jets  into  uniform  size  droplets,  resulting  in  a  monosized 

spray.  Precise  control  of  fuel  flow  rate  and  droplet  size 
is  obtained.  Preliminary  combustion  studies  demonstrate 
the  applicability  of  this  spray  generator  for  parametric 

studies. 


INTRODUCTION 

Methods  to  generate  a  controlled  spray  have  been  of 
major  concern  to  heterogeneous  combustion  research  for 
years.  Such  sprays  would  be  very  useful  in  that  the  effect 
of  important  parametars.  such  as  tht  droplet  size  and 
overall  fuel-oxidizer  ratio,  on  the  combustion  could  be 
checked.  Different  methods  have  bean  used  to  generate  a 
monodisptrse  spray.  Burgoyne  and  Cohen  (1,2)  used  a 
spray  generator  based  on  the  principle  of  Sincleir-LaMer  (3) 
and  LaMar  et  al.  (4)  aerosol  generator  to  study  the  effect 
of  droplet  size  on  the  burning  velocity  in  sprays.  The 
monosize  spray  was  obtained  by  condensation  of  saturated 
fuel  vapor.  The  simultaneous  control  of  the  droplet  size 
and  the  amount  of  condensed  fuel  is  extremely  difficult. 
Therefore,  air  fuel  ratio  in  Burgoyne  and  Cohen's 
experiments  was  not  a  constant  over  the  range  of  droplet 
sizes  studied.  Only  droplet  sizes  as  large  as  SO^m  can  be 
produced;  for  larger  droplet  sizes,  this  is  not  an  attractive 
method. 

Elaetrostatic  atomization  was  dtvtloped  by  Vonnegut 
and  Naubaur  (6)  to  produce  uniform  size  droplets  with 
diameter  of  l^m  or  less.  Rotating  disk  spray  generator 
has  baen  used  for  combustion  studies  by  many  researchers 
(6,7)  where  the  droplet  size  and  number  density  are 


NOMENCLATURE: 


D 

m 

diameter 

f 

m 

frequency 

n 

m 

number  of  orifices 

u 

■ 

velocity 

OD 

■ 

outside  diameter 

V 

m 

volume 

• 

density 

0 

■ 

surface  tension 

\ 

• 

wave  length 

4 

• 

sir/fuel  ratio 

SUBSCRTT* 
d  »  droplet 

i  •  J* 

I  •  liquid 

min  ■  minimum 

o  *  oxidizer 

opt  ■  optimum 

t  •  tube 


controlled  by  adjusting  the  rotating  speed  of  the  disk  and 
the  fuel  feed  rate,  respectively.  The  fuel  is  atomized  into 
droplets  in  all  directions  end  the  control  Of  droplet  flow 
without  loss  of  monodispersity  is  difficult.  Also  the  work 
of  Friedman  at  al.  (8)  indicated  that  wnen  these  sprays  are 
loaded  heavily,  the  uniformity  of  the  spray  diminishes. 

Vibrating  capillary  technique  which  was  developed  by 
Dimmoek  (9.10)  (10-300 ?  drops).  Margarvey  and  Taylor  (11) 
(S00^m-20mm).  and  Mason  et  al  (12.13)  (30-1000^)  to 
produce  a  stream  of  monosized  droplets,  was  modified  by 
Dabora  (14)  to  produce  a  monodispersad  spray.  He  used 
ten  vibrating  capillary  needles  fitted  in  a  plate  to  produce 
ten  streams  of  monosize  droplets.  Many  streams  are 
needed  to  represent  a  spray. 

The  impulsed  multi-orifice  spray  generator  developed 
here  suits  parametric  studies  of  spray  combustion.  This 
spray  generator  can  be  operated  with  two  different  modes. 
When  the  fuel  chamber  is  not  pressurized,  an  array  of 
droplets  are  ejected  from  the  orifices  on  the  plate  with 
each  pulse  of  the  piezoelectric  transducer.  The  axial 
spacing  of  the  droplet  streams  can  be  varied  directly  with 
the  frequency  of  the  pulse,  but  e  spectrum  of  droplet  sizes 
is  produced.  The  other  operating  method  of  the  same 
device  is  to  pressurize  the  fuel  chamber  in  which  droplets 
with  uniform  size  are  produced  with  the  optimum 
frequencies  of  Rayleigh  jet  breakup  by  the  piezoelectric 
oscillation. 

The  multi-orifice  plate  spray  generator  offers  the 
advantage  of  flexibility  in  parametric  studies  of  spray 
combustion.  In  this  paper  we  will  describe  the  details  of 
this  generator  and  its  operating  behavior.  Also  the  design 
of  en  apparatus  for  the  combustion  study  of  laminar 
monodispersad  droplet  flow  with  the  preliminary 
combustion  results  presented. 

MULTI-  ORIFICE  IMPULSED  SPRAY  GENERATOR 

The  principle  of  the  multi-orifice  impulsed  spray 
generator  design  is  to  utilize  the  pressure  pulse  from  a 
piezoelectric  transducer  (Model  104-95NS,  Piezoelectric 
Product  Inc.)  to  control  the  ejection  of  the  droplets.  A 
large  circular  piezoelectric  transducer  is  used  to  pressurize 
the  fuel  chamber  and  to  eject  the  droplets  from  a  plate 
with  many  orifices  on  it. 

The  generator  which  has  beeh  designed  is  composed 
of  a  body  (see  Fig.  1)  with  a  multi-orifice  plate  mounted  at 
the  top  side  and  a  piezoelectric  transducer  at  the  bottom 
aide.  The  interior  of  the  generator,  which  is  the  fuel 
chamber,  is  machined  smoothly  without  sharp  edges, 
because  sharp  corners  may  induce  cavities  during  pulsation. 
Since  the  displacement  of  the  piezoelectric  is  very  small, 
the  existence  of  a  small  cavity  would  absorb  the  pressure 
pulse  of  the  piezoelectric  transducer. 

Since  the  contraction  and  expansion  of  the 
piezoelectric  crystal  under  .electric  pulse  is  very  small, 
bilaminar  plate  arrangement  (15)  is  used  to  amplify  the 
displacement.  This  bilaminar  plate  is  made  by  epoxying 
the  piezoelectric  crystal  on  a  larger  brass  plate  of  4  cm 
diameter.  When  the  axially  polarized  piezoelectric  plate  is 
exposed  to  an  electric  field  in  the  direction  of  the 
polerization,  there  appears  an  axial  expansion  and  radial 
contraction  of  the  piezoelectric  plate.  The  axial  expansion 
is  small,  however,  the  radial  contraction  gives  a  concentric  1 
moment  that  bends  the  composite  plate  toward  the  fuel 
chamber  as  an  applied  pressure  impulse. 
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Figure  1  Pressurized  flow  operation.  Uniform  size 

droplet  streams  are  produced. 

Two  types  of  operations  of  this  impulsed  spray 
generator  can  be  performed: 


When  the  fuel  in  the  chamber  is  not  pressurized,  the 
droplets  can  be  ejected  from  all  the  orifices  at  any  applied 
frequency  of  the  piezoelectric  transducer.  When  a  voltage 
pulse  is  applied  to  the  piezoelectric  plate  •  pressure  pulse 
is  generated  in  the  container.  As  e  result,  en  array  of 
droplets  is  ejected.  When  the  pulse  ends,  the  plate  returns 
to  its  original  position  end  the  orifices  on  the  generator 
are  refilled  with  fuel  by  capillary  affect.  Pressure  inside 
the  chamber  is  very  close  to  the  atmospheric  pressure. 
The  initial  design  parameters,  such  as  the  orifice  length 
and  shape,  pulse  pressure,  and  pulse  duration  developed  by 
Beasly  (16)  for  a  tingle  orifice  was  modified  for  multi¬ 
orifice  conditions  in  (17). 

It  is  very  important  to  maintain  a  minimum  thickness 
of  the  orifice  plate  to  prevent  air  bubbles  from  being 
tucked  into  the  fuel  chamber.  This  problem  usually  occurs 
when  the  miniscus  retreats  back  into  the  chamber  after  the 
droplet  is  ejected.  A  thick  plate  (1.2  cm)  and  a  thin  plate 
(0.6  mm)  have  been  tested  and  the  thick  plate  appeared  to 
be  the  better  one  in  eliminating  the  etr  bubble  in  the 
chamber.  On  the  other  hand,  the  thicker  the  plate,  the 
more  the  viscous  force  will  attenuate  the  pressure  pulse. 
This  problem  can  ba  overcome  by  rfiaking  the  orifice 
diameter  varying  fs  a  converging  nozzle  viewing  from  the 
fuel  chamber. 

An  optimal  distance  between  the  piezoelectric  plate 
to  the  orifice  plate  is  required.  The  constraint  on  the 
piezoelectric  operation  is  the  maximum  applied  voltage. 
Therefore  maximum  pressure  and  displacement  also  occur 
for  aach  piezoalectric  plata.  In  ordtr  to  offer  the  required 
magnitude  of  the  pressure  pulse  to  the  orifices,  a  short 
distance  between  the  piezoelectric  plate  and  the  multi¬ 
orifice  plate  is  desirable.  However,  very  short  distance 
would  generate  an  array  of  droplets  with  non-uniform 
sizes.  The  droplets  near  the  center  of  the  plate  would  be 
ejected  faster  due  to  the  uneven  buckling  of  the 
piezoelectric  plate. 

In  this  type  of  operation,  the  droplet  generation 


» 
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frequency  can  t>e  varied  at  des  ired.  However,  the 
generated  droplets  usually  show  a  wide  spectrum  of  sizes. 
This  is  because  the  orifices  at  different  positions  on  the 
plate  will  experience  different  magnitudes  of  pressure 
pulses.  In  practice,  if  there  exists  a  liquid  film  on  the  top 
of  the  orifice  plate,  the  generated  droplet  will  have  to  tear 
away  from  the  liouid  film  and  will  have  a  different  size 
from  that  generated  from  a  hole  on  a  dry  plate.  This 
liquid  film  exists  usually  due  to  the  minor  leakage  of 
liquid  from  the  slight  pressure  differences  across  the 
orifice  plate.  Additionally  any  imperfections  of  the  orifice 
hole  will  cause  the  disturbance  of  the  pressure  wave  in  the 
hole  and  the  generation  of  satellite  droplets  during  the 
injection. 

Pressurized  Flow  Operation 


The  other  type  of  operation  of  the  same  device  is  to 
pressurize  the  fuel  chamber  so  that  liquid  jets  flow 
continuously.  With  a  specific  pulse  frequency  applied  at  a 
jet  velocity,  the  liquid  jet  breaks  up  into  uniform  size 
droplets  (Fig.  1.).  The  phenomenon  of  uniform  liquid  jet 
break-up  was  first  observed  experimentally  by  Savart  (18) 
and  studied  theoretically  by  Plateau  (19)  and  Rayleigh 
(20,21).  From  the  consideration  of  surface  energy,  Plateau 
derived  the  minimum  wavelength  X  (distance  between 
unstable  disturbance)  on  an  infinitely  long  liquid  cylinder. 


where  D  is  the  liquid  jet  diameter.  For  an  inviscid. 
incompressible,  cylindrical  liquid  jet  sprayed  into  a  vacuum 
Rayleigh  derived  the  optimum  wavelength  for  a  most 
unstabla  disturbance. 

•  « 08  D,  (2) 
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Schneider  end  Hendricks  (22)  experimentally 
that  uniform  droplets  could  ba  produced  by 
varying  X  within  the  following  range: 

X5  DJ  <  X  <  7  Dj  (3) 

Wabar  (23)  extended  the  analysis  to  include  the 
effects  of  the  liquid  viscosity  and  the  velocity  of  the  jet. 
but  Crane  at  al  (24)  showed  that  if  the  velocity  is  low 
(less  than  10  m/()  these  effects  can  be  neglected. 

Considering  the  conservation  of  energy  Lindblad  end 
Schneider  (25)  derived  e  minimum  liquid  velocity  needed  to 
form  s  liquid  jet  from  a  cspillery  tube. 


ui  ■  ^~T^U2  (4) 

where  9  is  the  surface  tension  and  p  is  tha  dansity 
of  the  liquid.  Dsbora  (14)  experimentally  showad  that 
actual  minimum  valocities  art  lowar  by  25*35%  than  this 
relation  indicated. 


The  diameter  of  the  liquid  droplet  particles  from  the 
impulsed  orifice  generator  can  ba  calculated  from  jat 
velocity  and  the  disturbance  frequency  using  simple  mass 
conservation  aquation.  That  is 

Od  -  D  *  u /f],/3.  (5) 

This  kind  of  operstion  is  vary  similar  to  that  of  the 
Berglund  81  Liu  Aerosol  Generator  (26);  however,  the 
present  device  generates  droplets  from  many  orifices  so 
that  a  spray  with  very  dense  droplets  end  little  turbulence 
is  produced.  The  number  of  active  orifices  on  the  plate 
can  also  ba  varied  by  plugging  tha  orifices  or  by  using  a 


different  plate  with  different  numbers  and  patterns  of 
orifices.  Therefore,  to  generate  a  spray  with  a  fixed 
droplet  size  from  a  given  orifice  plate  the  u  If  is  changed 
(see  equation  (5))  to  give  different  fuel  flowJ  rates.  When 
the  operation  beyond  the  limits  of  equation  (3)  is  desired, 
the  plate  is  substituted  with  another  plate  with  different 
orifices. 


The  number  of  operating  orifices  can  be  related  to 
the  air/fuel  ratio,  f  for  any  particular  droplet  size  and  air 
flow  velocity  from  the  consideration  of  mass  conservation. 


3  °,2  u„ 
n  •  -  — »  — 2  • 

2  D?  f 


(6) 


-6  ■ 

As  an  axampla,  for  a  stoichimetric  hexane  air  mixture 
with  2  m/s  air  valocity  and  dropfat  of  285^m  at  3000HZ, 
12  orifices  era  needed  on  a  plate  with  5  cm  diameter. 


COMBUSTION  FACILITY 

In  order  to  study  the  combustion  behavior  of  a  spray 
generated  by  the  present  device,  a  special  experimental  set 
up  is  fabricated  as  shown  in  Figure  2.  Generated  droplets 
flow  vertically  upward  through  several  pieces  of  tubes 
which  have  a  diameter  of  9  cm  and  length  of  4  cm  each. 
A  uniform  suction  at  the  downstream  moves  the  droplet 
flow  vertically  upward.  Air  is  sucked  from  the  quiescent 
atmosphere  through  the  spacings  of  the  tubes.  As  a 
result,  the  droplets  will  not  collide  with  the  well. 


Air  is  also  ejected  into  the  droplet  flow  through  an 
array  of  tubes  situated  in  between  the  rows  of  the  orifices 
on  the  plate.  The  air  tubas  are  0.3 1  cm  OD  with  aligned 
holes  drilled  on  them.  The  hole  positions  are  adjusted  to 
give  uniform  flow,  the  pressurized  air  flow  is  metered 
before  fed  into  the  tubes.  With  the  air  jet  flow,  the 
stream  of  droplets  may  be  slightly  dispersed  to  prevent 
coalescence  of  the  droplets. 

The  spray  is  ignited  with  a  hydrogen  pilot  flame. 
Very  small  amounts  of  hydrogen  flow  is  ejected  from  the 
fine  holes  on  the  top  side  of  a  horizontal  tube  of  1.6  mm 
O.D.  After  ignition,  the  hydrogen  flow  can  be  terminated 
then  the  tube  also  acts  as  a  flame  holder.  A  steady  V 
shaped  spray  flame  can  be  sustained  on  this  tube. 
Knowing  the  gas  and  liquid  flow  velocities,  and  the  stream 
lines  flame  speed  can  be  estimated  from  the  angle  of  the 
V  flame. 


RESULTS  AND  DISCUSSION 

The  characteristics  of  the  multi-orifice  impulsed  spray 
generator  are  studied  for  two  different  multi-orifice  plates. 
In  making  the  60^m  orifice  plate,  •  special  carbide  drill 
has  been  used.  Table  (1)  indicates  the  dimension  and  the 
number  of  orifices  as  well  es  the  material  of  each  plate. 
The  estimated  size  variations  of  the  orifice  is  within  5  p m. 
Also  the  ranges  of  the  monosize  droplets  at  pressurized 
flow  operation  are  also  indicated  in  the  Table. 
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Figure  2 


Overall  Setup  for  the  Combustion  Studies. 


Table  1 


Multiorifice  Plate 
Orifice  Diameter 

#1 

42 

Exit 

60  ^m 

160 

Inlet 

0.6  mm 

2  mm 

Number  of  Orifices 

221 

39 

Orifice  Spacing 

3x3 

4x10 

Plate  Thickness 

0.6  mm 

1.2  cm 

Materiel 

Stainless  steal 

Brass 

Orifice  Nozzle 

Connicai 

Stapwise 

Shape 

Cylinders 

Plate  Diameter 

6.9  cm 

6.9  cm 

Droplet  Diem. 

Range  Obtainable 

104-131 

260-330 

Figure  3  shows  the  spectrum  of  droplet  size  as  the 
percent  of  volume  for  I50^m  orifice  plate  under  impulsed 
jet  operating  conditions.  A  microscope  with  magnification 
up  to  112  times  is  used  for  the  photographing  of  the 
droplet  sizes.  The  strobe  light  is  triggered  by  the  shutter 
of  a  Polaroid  camera  which  is  attached  to  the  microscope. 
Droplet  size  spectrum  in  Figure  3  is  found  from  the 
measurement  over  1000  droplets  from  photographs  that 
were  magnified  14  times.  Although  a  wide  spectrum  is 
observed,  it  is  clear  that  there  is  a  peak  between  220  to 
300/m.  The  smaller  peak  near  I50^m  is  possibly 
originated  from  the  satellite  droplets.  The  spray  generated 
at  this  operation  condition  is  suitable  for  studying  poly* 
dispersed  spray  combustion  with  the  droplet  spacing 
adjustable  arbitrarily. 

Figure  4  shows  the  size  spectrums  of  the  droplets  at 
pressurized  flow  operating  condition  for  I50^m  orifice 
plate  at  two  different  flows  and  frequencies.  Monosize 
droplet  spray  with  about  1  10^m  deviation  from  the  mean 
droplet  size  can  be  obtained.  Fig.  <4-6).  This  deviation  is 
possibly  because  the  orifices  are  not  exactly  the  same  in 
size.  Also  the  u  is  possibly  slightly  different  for  the 
orifices  at  the  sidis  end  at  the  middle  of  the  plate.  This 


is  observed  by  the  maximum  height  droplets  reach  when 
ejected  vertically  up  under  stagnant  ambient.  It  is  seen 
that  the  maximum  height  reached  by  droplets  ejected  from 
Sice  orifices  is  »lightly  smaller  than  the  ones  in  the  center 
of  the  plate.  Larger  sine  droplets  will  be  produced  near 
the  middle  of  the  plate  if  the  averaged  pressure  is  higher 
there.  The  spectrum  in  Figure  4  is  obtained  by 
photography  measurement.  The  typical  prieute  of  the 
droplet  streams  taken  close  to  the  plate  is  shown  in  Figure 
5. 


100  140  ISO  220  260  300  340  360  420 


Figure  3  Droplet  size  spectrum  by  photography 

measurement  technique,  for  impulse  jet 
ejection  operation  at  1000H2  frequency. 
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Figure  4  Droplet  size  spectrum  by  photography 

measurement  technique  for  pressurized 
flow  operation.  (A.  3700HZ  and 

40cc/min(  B.  at  2800HZ  and  40cc/min 
Hexane  Flow  Rate.) 


When  spray  is  generated  at  pressurized  flow 
condition,  the  required  fuel  flow  rate  and  the  impulse 
frequency  to  give  a  desired  droplet  size  can  t>e  evaluated 
prior  to  the  experiment  with  good  accuracy.  Figure  6 
shows  the  operating  range  of  the  I50ym  mu'ti-orif  ice 
plate.  The  obtainable  range  of  drop'et  size  is  between 
260,/m  and  330ym.  This  can  be  evaluated  from  equations 
(3)  and  {5)  assuming  3k  *  u  /f  where  the  X  is  variable  in  a 
range  as  described  in  equation  (33.  The  curves  of  droplet 
size  versus  frequency  are  drawn  using  equation  (5)  with 
three  different  jet  velocities  of  2.1.  2.6  and  3.3  ml  as  the 
solid  line  in  Figure  (6).  These  jet  velocities  are  evaluated 
from  the  measured  mass  flow  of  fuel  through  the 
generator.  The  hexane  flow  rate  can  be  set  to  produce  the 
required  jet  velocity.  Then  droplet  sizes  at  different 
frequencies  are  obtainable.  The  agreement  between  the 
experimental  points  and  the  curves  predicted  from  equation 
(5)  is  reasonable. 


Streams  of  uniform  size  droplets  near  the 
orifice  plate. 


In  order  to  demonstrate  the  combustion,  a  monosce 
spray  of  285ym  mean  droplet  diameter  generated  at  30:0 
HZ  frequency.  2.1  m/  jet  velocity  from  19  oc^rat.ng 

orifices  is  considered.  Air  tubes  are  adjusted  to  a  ve  a 
steady  upward  homogeneous  flow  of  droplets  /j  the 
ign;;ion  point.  Gas  velocity  profiles  are  mtasured 

separately  under  cold  condition  with  no  droplet  flow  ysing 
a  hot  film  anemometer. 

Pictures  of  the  spray  are  taken  near  the  ignition  pomt 
using  strobe  light  and  the  microscope.  Examining  the 
photographs,  it  is  noticed  that  the  spray  has  always  kept 
its  monosize  distribution.  Pictures  of  the  sustained  V- 
shaped  flame  on  the  horizontal  flame  holder  are  shown  ir 
Figure  7(a).  Figure  7(a)  shows  the  flame  of  isohexane 
spray  with  initial  droplet  size  of  about  285  ,m.  By 
studying  this  steady  flame,  detailed  behavior  of  the  flame 
propagation  in  the  vicinity  of  the  flame  front  can  be 

observed  with  the  microscope.  A  spark  source  (strobe 

light)  with  3  ft s  duration  is  used  to  illuminate  droplets. 
Figure  7(b)  shows  the  flame  of  isohexane  spray  with  imtiai 
droplet  size  of  150  t  10  ,m.  It  is  shown  that  all  droplets 
vaporize  completely  ahead  of  the  flame  and  a  blue  flame 
is  obtained  And  Figure  7(c)  shows  the  flame  of  Toluene 
spray  with  initial  droplet  size  of  285  ym.  No  blue  flame 
is  observed. 

Through  various  combustion  studies,  it  is  observed 
that  the  large  droplets  of  285ym  would  not  completely 
vaporize  in  the  preheat  zone.  When  the  droplet  number 
density  is  low,  individual  droplets  burn  in  blue  envelope 
flames.  As  the  number  density  is  increased,  the  flame 
shape  and  structure  changes  to  a  yellow-  flame  with  a 
rough  flame  front.  The  details  of  the  combustion  results 
will  be  presented  in  a  future  paper. 


spark  source.  Spark  light  show-s  droplets 
approaching  the  flame  front. 


Figure  7b  Isohexane  spray  flame,  with  initial  droplet 

size  of  130  fitr\ .  The  blue  flame  shows 
ell  the  droplets  are  vaporized  ahead  of 
the  flame. 


Figure  7c  Toluene  spray  flame  with  initial  droplet 

size  of  285  No  blue  flame.  Flame 

front  is  not  smooth. 


CONCLUSION 

A  multi-orifice  impulsed  spray  generator  with  two 
different  modes  of  operation  was  developed.  The  impulse 
jet  operation  mode  covers  a  wide  range  of  droplet 
frequency  and  produces  arrays  of  droplets  with  controllable 
Spacing,  but  with  a  wide  s,:e  spectrum.  Pressurized  flow 
operation  applies  disturbance  frequency  to  break  up  the 
liquid  jets  into  uniform  size  droplets,  resulting  in  a 
monosized  spray.  Precise  control  of  fuel  flow  rate  and 
droplet  size  is  obtained.  Preliminary  combustion  studies 
demonstrated  the  applicability  of  this  spray  generator  for 
parametric  studies. 


REFERENCES 

1.  Burgoyne,  J.H.  and  Cohen.  L.,  "Effect  of  Orop 
Size  on  Flame  Propagation  in  Liquid  Aerosols." 
Prpc.  Roy.  Soc.,  1954.  Vol.  A225,  No.  1162.  pp. 
375-392. 

2.  Burgoyne.  J.H.  and  Cohen.  L,  "Production  of 
Monodispersed  Aerosols  of  Large  Drops."  J. 
Colloid  Scu.  1953.  8(3).  p.  364. 

3.  Sinclair.  D.  LaMer,  VX,  "Light  Scattering  as  a 
Measurement  of  Particle  size  in  Aerosols." 
Chem.  Rev.,  1949,  44.  p.  245. 

4.  LaMer.  VX,  Inn.  E.C.Y.  and  Wilson,  I.B.,  "The 
Method  of  Forming,  Detecting  and  Measuring  the 
Size  and  Concentration  of  Liquid  Aerosols  in 
the  Size  Range  0.01-0.25  Microns  Diameter,"  J. 
Colloid  Sci.,  1950.  S,  pp.  471-497. 

5.  Vonnegut.  B.  and  Neubaur,  R..  "Production  of 
Monodispersed  Liquid  Particles  by  Eiectral 
Atomization."  J.  Colloid  Sci..  1952.  7(6),  p.  616. 

6.  Hayashi,  S.,  Obtani,  T.,  Iinuma.  K..  and  Kumagai. 
S.,  "Limiting  Factors  of  Flame  Propagation  in 
Low-Volatility  Fuel  Clouds,"  Eighteenth 
Symposium  (Int.)  on  Combustion.  The 
Combustion  Institute  1981,  pp.  361-367. 

7.  Bolt,  J.A.  and  Boyle,  T.A.,  "The  Combustion  of 
Liquid  Fuel  Spray,"  Trans.  ASME  78.  (1956), 
609. 

8.  Friedman,  S.J.,  Gluckert,  F.A.  and  Marshall.  W.R., 
Jr.,  Chemical  Engineering  Progress.  Vol.  48, 
April,  1952,  pp.  181-191. 

9.  Dimmock,  N.A.,  "Production  of  Uniform 
Droplets,"  Nature  (Engl),  1950.  Vol.  166  pp.  $86- 
687. 

10.  Dimmock,  N.A.,  "The  Controlled  Production  of 
Streams  of  Identical  Droplets,"  Nat.  Gas  Tub. 
Est.  (Engl),  Memo  M  115,  1961. 

11.  Magarvey,  R.H.  and  Taylor.  B.W..  "Apparatus  for 
the  Production  of  large  Water  Drops/’  The  Rev, 
of  Sci.  Inst..  Nov.  1956,  Vol.  27,  No.  11.  pp. 

944-47. 

12.  Mason,  B.J.,  Jayoratne,  P.W.,  and  Woods,  J.D., 
"An  Improved  Vibrating  Capillary  Device  for 
Producing  Uniform  Water  Droplets  of  15  to 
500^m  Radius,"  J.  Sci.  Inst..  1963.  Vol.  40.  pp. 
247-49. 

13.  Mason,  B.J.  and  Brownscombe,  J.L.,  "Production 
of  Uniform  Size  Drops  at  Controllable  Frequency 
and  Spacing  from  a  Vibrating  Capillary.”  J.  Set. 
Inst..  May  1964,  Vol.  41,  pp.  258-259. 

14.  Dabora,  EX,  "Production  of  Monodispersed 
Sprays."  The  Review  of  Scientific  instruments. 
Vol.  38,  No.  4,  1967. 

15.  Ashton,  J.E.  and  Whitney  J.M.,  "Theory  of 
Laminated  Plates,"  Stanford,  Conn.:  Techncmic, 
1970,  pp.  1-29. 


16.  Beasly.  J.D.,  "Model  for  Fluid  Ejection  and  Refill 
in  Air  Impulse  Drive  Jet,"  Photographic  Science 
and  Eng..  Vol.  21.  No.  2  (March/April,  1977). 

17.  Athgrizzadah.  M.  Yao.  5.C,  and  Sirigneno,  W.A., 
"Progress  on  Impulse  Jet  Droplet  Generator  for 
Laminar  Spray  Combustion  Experiments,"  Eastern 
Section/Combustion  Institute.  Princeton,  New 
Jersey,  p.  44,  (1980). 

16.  Savart.  Ann.  Chimie  53,  337  (1833). 

19.  Plateau.  "Statique  Experimentale  at  Tehorique 
des  Liquids  aaumis  aux  seules  Forces 
Moleculaires."  (1873)  Ref.  in  "Theory  of  Sound." 
by  J.W.S.  Rayleigh.  2nd  ed..  Vol.  11.  363  (1878). 
Reprinted  by  Dover  Publi.,  N.Y.  1945. 

20.  Rayleigh,  Lord,  Proc.  Roy.  Soc.  29.  71.  (1879). 

21.  Rayleigh,  Lord.  Proc.  London  Math.  Soc..  10.4, 
(1878). 

22.  Schneider.  J.M.,  Hendricks.  C.D-  Rev.  Sci. 
Instrum.  35  (10).  1349-50  (1964). 

23.  Weber,  C.,  2.  Angew.,  Math  and  Mech.  11,  136. 
(1931). 

24.  Crane,  L..  Birch,  S.  and  McCormack,  P.D..  Brit.  J. 
Appl.  Phys.  15.  743,  (1964). 

25.  Lindbiad,  N.R.  and  Schneider,  J.M..  "Production 
of  Uniform-sized  Liquid  Droplets,"  J.  Sci.  tnstr. 
42,  635,  (1965). 

26.  Berglund.  R.N.  and  Liu.  Y.H..  Environmental  Sci. 
6t  Tech.,  Vol.  7,  pp.  107-153,  (1973). 


CFFECT  CF  THERMAL  RADIATION  Of!  THE 
DROPLET  PRE-EVAPORATION  If,1  SPRAY  COMBUSTION 

Shi-chune  Yao  flassar  Ashgri zzadeh 

Associate  Professor  Graduate  Student 


Department  of  Mechanical  Engineering 
Carnegi e-Mel  Ion  University 
Pittsburgh,  PA  15Z13 


INTRODUCTION 

In  many  combustion  systems  the  flame  and  the  combustion  products  are  at  very 
high  temperatures  such  that  thermal  radiation  plays  a  significant  role  in  the  trans 
port  of  combustion  heat  [1].  Generally,  the  combustion  process  of  gases  is  not 
strongly  influenced  by  che  thermal  radiation  from  the  hot  zcr.e  because  the  abserpti 
vitics  of  gases  are  usually  small.  However,  the  combustion  process  of  two  phase 
mixtures  could  be  affected  by  radiation  significantly  because  of  the  high  absorp- 
tivities  of  the  participating  liquid  or  solid  phases. 

Different  from  heat  convection,  thermal  radiation  is  a  far-field  interaction. 


In  spray  combustion  the  droolets  could  be  heated  up  dnd  possibly  evaporated  by  the 
radiation  at  the  far  upstream  cf  the  flame.  Before  the  spray  reaches  the  flame, 
the  droplet  size .diminishes  ,  the  droplet  number  density  reduces,  and  the  fuel- 
vapor  to  air  ratio  increases  as  compared  with  the  equivalent  case  where  radiation 
does  not  exist.  Therefore,  the  subsequent  combustion  of  the  spray  is  greatly 
affected  because  the  thermal  radiation  changes  the  initial  condition  of  the  spray 
combustion.  In  extreme  cases,  the  droplets  may  be  pre-evaporated  completely  and 
the  combustion  of  the  spray  becomes  a  combustion  of  pre -mixed  gases. 

The  objective  of  this  pacer  is  to  study  the  radiative  effect  on  the  droplet 


flow  combustion  to  reveal  the  fundamental  behavior  of  th*s  interaction. 


radiative  transport  will  be  calculated  for  an  idealized  one  dimensional  droplet 
flow.  The  transfer  number  B  of  droplet  considering  the  thermal  radiation  is 
calculated.  Finally,  calculation  is  performed  for  a  sample  case  to  demonstrate 
quantitatively  the  radiative  effect  on  droplet  pre-evaporation. 


MODEL 

In  order  to  demonstrate  the  basic  principle  of  radiative  transport,  a  simple 
model  of  one-dimensional  droplet  flow  combustion  system  is  considered  with  its 
schematic,  shown  in  Figure  1. 


With  a  stationary  flame  in  space,  the  combustion  products  constitute  a  hot 
zone  with  a  finite  depth  and  at  uniform  temperature.  Monodisperse  droplets  are 
suspended  at  the  far  upstream.  Due  to  radiative  evaporation  the  droplet  diameter 
d  becomes  a  function  of  space  r.  The  hot  gases  perform  as  radiation  source  with 
self  absorption  while  the  droplets  at  upstream  of  flame  perform  as  radiative  sink. 

III.  RADIATIVE  TRANSPORT 

It  is  intended  to  find  the  radiative  heat  flux  to  a  droplet  of  diameter  d  at 
a  distance  h  upstream  of  the  flame.  The  droplet  is  spherical;  therefore,  the 
radiation  from  any  angle  will  see  it  as  a  round  disk  with  an  area  a-d^r.  The  total 
radiative  heat  flux  qxd  for  the  wave  length  A  from  all  the  differential  areas  j  to 
the  droplet  is  — 

At  Ut  -  f  e“)  dk  AJ  fJ*  r-U  01 

where  e^  is  the  black  body  hemispherical  emissive  power  and  (aiT2 ) jd  9e0- 

metric  mean  value  of  the  emittance  in  hot  region  and  transmittance  in  the  droplet 
flow  region  from  differential  area  j  to  droplet  d.  Considering  the  integration  of 
all  the  symmetric  rings  of  surface  element  j  the  equation  (1)  gives  the  incident 
heat  flux  to  the  droplet  at  the  wavelength  A 
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For  an  one-dimensional  system  the  0  varies  from  0  to  7  .  Integrating  the 
equation  (2)  it  becomes 

jkA  =  ^  [ez  t)-£,  OU.H  +  a,„  i  ) J  ek6j  va 


where  the  E2  is  the  second  exponential  integral  function. 

Not  all  the  radiation  incident  to  the  droplet  is  absorbed  by  the  droplet. 
Only  a  fraction,  which  equals  the  absorption  efficiency,  will  be  absorbed.  In  the 
above  analysis  it  has  been  assumed  that  forward  scattering  occurs  to  the  droplets 
and  the  radiation  property  of  the  two  phase  mixture  is  optically  thin.  For 
simplicity,  it  is  also  assumed  that  the  droplets  burn  instantaneously  when  they 
go  into  the  flame  and  the  burning  droplets  in  the  hot  region  do  not  participate 
in  radiative  energy  transport. 

In  many  practical  conditions  the  above  assumptions  are  satisfied.  For 
dielectric  fuel  droplets  with  diameter  larger  than  Bum  the  radiation  from  hot  gas 
at  a  temperature  higher  than  1500  K  will  have  predominantly  forward  scattering. 

The  optically  thin  approximation  requires  that  the  photon  mean  free  path  is  much 
larger  than  the  characteristic  dimension  of  the  system.  For  droplets  with  ICCum 
diameter  and  2mm  interspacing  the  radiation  properties  car.  be  considered  optically 
thin  with  respect  to  a  characteristic  dimension  less  than  lm.[2]. 

TRANSFER  NUMBER  FOR  RADIATIVE  EVAPORATING  DROPLET 

The  classical  transfer  numbers  of  evaporating  or  combusting  droplet  have 
been  well  established  for  conductive  or  convective  heat  transfer  conditions. 
However,  these  transfer  numbers  should  be  revised  for  cases  where  radiation  is 
significant. 

Since  the  gas  phase  does  not  participate  actively  in  the  radiative  transfer, 
the  transfer  number  of  a  droplet  can  be  rederived  by  simply  changing  the  boundary 
condition  of  the  vapor  energy  equation.  [3] 
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By  performing  integration,  the  transfer  number  considering  radiative  effects 
becomes  -  *  _  . 

8  =  -  Cr>CT„- Tc) 


(4) 


where 
species  equation 


is  the  vapor  velocity  at  the  droplet  surface  which  can  be  derived  from 


24  jr  =  C/*e„)/rs  (s) 

where  the  form  of  3^  is  not  changed  due  to  radiative.  It  is  interesting  tc  point 
out  that  the  transfer  number  By  is  therefore  indirectly  dependent  upon  the  droplet 
diameter  rs.  As  a  result,  the  diminishing  of  droplet  size  will  no  longer  follow 
the  dz  rule  exactly. 


AUXILLIARY  EQUATIONS  AND  SAMPLE  CALCULATIONS 

The  change  of  droplet  size  due  to  radiative  pre-evaporation  can  be  calculated 
when  the  B  value  is  known.  The  B  will  be  evaluated  from  the  B  value  which  is,  in 
turn,  dependent  upon  the  radiative  heat  flux  from  the  hct  zone.  In  general,  the 
radiative  flux  is  a  function  of  x  and  the  size-and-distribution  of  droplets  between 
the  location  x  and  the  flame  front.  Conservation  equations  can  be  written  for  the 
evaluation  of  this  information  without  substantial  difficulty;  however,  the  calcu¬ 
lation  will  be  performed  numerically. 

In  spite  of  the  practical  complications  of  the  auxilliary  equations,  sample 
hand  calculation  can  be  performed  to  reveal  the  essential  feature  of  this  problem 
with  some  approximations.'  For  example,  the  radiative  transport  can  be  calculated 
easily  if  a  gray  system  is  assumed.  For  optically  thin  two-phase-mixture  the 
local  heat  flux  will  net  be  sensitive  to  the  variations  of  droplet  size  and  distri¬ 
bution  in  the  space.  Therefore  the  d  and  n  can  be  fixed  as  constants  in  the 
approximated  method  of  solution. 

Typical  calculations  have  been  performed  for  stoichiometric  mixture  of  Cetane 
with  droplet  size  lOOym  at  latm.  The  flame  speed  is  assumed  to  be  40cm/sec  and 
the  length  H  of  the  hot  zone  is  selected  as  30cm.  The  calculated  transfer  number 
is  about  9.6  instead  of  5.0  which  is  the  typical  value  when  combustion  without  raoi 
ation  The  transfer  number  is  found  to  be  sensitive  to  the  thickness  and  tempera¬ 
ture  of  the  hot  zone  where  radiation  origin  from.  'As  a  result,  a  droplet  at 
6cm  upstream  of  flame  may  diminish  from  100pm  to  75ym  in  a  traveling  distance 
much  less  than  half  cm.  Of  course  at  other  cases  of  lower  flame  temperatures  the 
corresponding  distance  could  be  longer. 


NOMENCLATURE 


.a,jj  absorption  coefficient  of  the  combustion  product  at  hot  zone 

a„2  absorption  coefficient  of  the  two  phase  mixture 

By  transfer  number  due  to  thermal  effect 

transfer  number  due  to  mass  transfer  effect 
d  droplet  diameter 

h  location  of  droplet,  see  Fig.  1 

H  thickness  of  hot  zone,  see  Fig.  1 

L  latent  heat  of  evaporation  of  fuel 

Pd  radiation  heat  flux  to  the  droplet,  all  the  wavelength 
rs  instantaneous  droplet  radius 

B  evaporation  parameter  in  the  equation  for  d2  rate 
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1.  Introduction 

Radiation  effect  on  the  droplet  pre-evaporation  in  spray  combustion  has  recently  become 
a  subject  of  concern,  as  the  theoretical  models  for  droplet  spray  combustion  beeing  developed 
in  a  more  detailed  fasion.  Radiation  effect  on  evaporation  of  a  single  droplet  and  droplet 
sprays  has  usually  been  neglected,  because  of  its  small  influence  in  droplet  evaporation.  But 
most  experiments  in  this  area  were  done  on  pTotype  laboratory  scaled  systems.  Where  small 
scale  combustion  systems  were  being  studied.  So  analytical  results  without  consideration  of 
radiation  are  close  to  experimental  data.  Although  radiation  is  a  dominant  heat  source  in  large 
scale  practical  systems  and  should  be  included  in  the  analytical  models. 

In  this  paper  the  authors  have  followed  their  previous  work  in  Droplet  pre-evaporation 
due  to  radiation  Cl  J.  Previously  a  model  of  infinitly  long  flame  was  considered  for  its 

simplisity  of  radiation  calculations.  In  this  paper  a  more  practical  model  of  cylinderical 
combustion  chamber  is  considered.  Also  wall  effects  which  was  neglected  in  the  previous 

work,  is  included  in  present  paper.  Spaldings  approach  of  evaporation  of  a  single  fuel  droplet 
in  a  nonconvective  atmosphere  of  a  given  temperature  and  pressure,  which  was  used  in  previous 
paper [2],  is  substituted  with  a  better  model  of  spacially  uniform  but  temporally  varying 
droplet  temperature  developed  by  Law [3]. 

2.  Droplet  Evaporation  In  A  Radiative  Environment 

By  assuming  Fick’s  Law,  Fourier’s  Law  and  unitary  Lewis-Semenov  number  we  can  solve 
for  droplet  evaporation  [3.4].  A  model  of  spacially  uniform  but  temporally  varying  droplet 
temperature  in  a  non-convective  environment  is  solved.  Radiation  heat  is  assumed  to  be  the 
only  heat  source,  to  show  only  the  radiation  effect  on  the  droplet  evaporation.  For  a  typical 
n-Decane  fuel  droplet  of  lOO^m  diameter,  in  a  non-convective  surrounding  the  assumption  of 
uniform  droplet  temperature  is  good  as  the  droplet  life  time{O(10)msec)  is  much  larger  than 
droplet  heat  up  time{O(10"')msecl.  Also  for  simplicity  we  assumed  all  the  radiation  heat 

deposites  at  the  surface  of  the  droplet.  This  assumption  can  be  removed  at  the  expense  of 

incurring  additional  complexity  of  absorption  of  radiation  heat  at  different  depths  within  the 
droplets [5-7],  The  governing  equations  of  a  evaporating  droplet  with  above  assumptions,  is 
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well  established,  neglecting  the  radiation  heat[3l-  So  considering  radiation,  modification  of  the 
analysis  would  be  only  at  the  boundry  condition  of  the  droplet.  A  heat  balance  at  the  droplet 
surface  is: 

k(^L)=mL+i-r/>c(— )-q  (1) 

dr  '  »  v  3  i  '  dt  ,  R 


A  new  transfer  number  is  defined  as  follow: 


C(T  -  T) 

B'  =  _ E _ 5 - 


where  is  defined  as:  q^*-  (r^  C)  (dT/dt)#  is  the  heat  required  to  raise  the  droplet 


temperature. 


The  evaporation  rate  m  can  also  be  found  from  species  equation  as: 
p  D 

m  =  -i — s  ln(l+B  )  (3) 

r  m 

ft 

Numerical  calculations  for  B  values  show  that  B  in  despite  of  what  was  thought  previously, 
never  blows  up.  As  [L^  -  (qR  -  q^l/m]  does  not  go  to  zero  for  all  practical  values  of  qR. 
This  can  be  explained  by  the  direct  relation  between  radiative  heat  and  vaporization  rate.  So 
the  more  radiation  the  higher  the  vaporization  rate  m  .  So  B  value  defined  as  above  can  still 
be  a  valuable  parameter  in  problems  including  radiation  heat  transfer. 


3.  Radiation  Heat  Transfer 

In  order  to  find  radiative  heat,  we  assume  forward  scattering  to  the  droplets  and  the 
radiation  property  of  the  two  phase  mixture  is  optically  thin.  For  simplicity  it  is  also  assumed 
that  the  droplets  burn  instantenously  when  they  go  into  the  flame  and  the  bumning  droplets  in 
the  hot  region  do  not  participate  in  the  radiative  energy  transport.  A  better  explanation  for 
these  assumption  is  given  by  Yao  and  Ashgrizzadehtl]  .  The  radiation  comming  from  Hot 
burning  gas  region  with  considering  the  wall  effect  of  that  region,  to  a  droplet  located  in  the 
centerline  of  the  cylinderical  chamber  is  given  by [8]: 

dQ.  *  A  dq.  =  ^  .  (dq.  AF  r.  e  d.  AF  «.  )r. 

'^Xi.k  k  ^Xi.k  MXo,j  j  X.j-k(l)  |0>bt  X  j  j-k  Xj-k(t)  X,j-k(2> 
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Where  dq.  =  «  <r  T4  is  the  radiation  from  the  wall.  And  the  second  part  of  the  right 

A  w  w 

hand  side  S'  radiation  due  to  the  hot  gasses.  Both  terms  are  attenuated  by  the  factor  while 
passing  through  cold  zone  due  to  existance  of  absorbing  droplets  in  the  radiation  path. 

F  =(Cos/S Cos/?  )/r(S  +S  )2  (5) 

j-fc  j  k  12 

is  the  view  factor.  As  droplets  are  spherical,  the  radiation  from  any  angle  will  see  it  as  a 
round  disk  with  an  area  *d2/4  and  B  =0.0  .  Also: 

r  k 

a .  =1-  exp(-a.  S )  (6) 

X,j-ku>  ^  Xt  i 

r .  =1-  exp(-a.  S )  (7) 

X.j-k(2)  ^  X2  2 

r .  =1-  exp(-a.  S )  (8) 

X.rkw  r  Xi  i 

Integration  of  equation(l),  using  equations(2-4)  results  in  the  following  equation  for  radiation 
arriving  to  a  droplet  located  at  a  distance  h  from  the  flame. 


*>x  -  2’r{  E,(*x,6)-  *x,h>*  V^x,0*  •xx"'’1- 


•2(aXj  h2+RJ  -  (1-Iw)  3  h  expHc^h/S^  1 


where:  n=  1+ 


(h^D") 


and  radiation  from  the  wall  is  assumed  to  be: 


dq.  *  t  o  T 

A  wvrw 


The  integration  limits  were  from  S}*h  to  h2+  R2  .  The  last  part  of  equation(5)  is 
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integrated  numerically  using  Simpson’s  rule.  Also  it  was  noticed  that  radiation  from  the  walls 
are  very  important,  and  if  combustor  has  long  wall  a  greate  proportion  of  radiation  heat  is 
supplied  by  the  wall.  Fig  2  shows  the  radiation  flux  to  droplets  with  diameter  of  100 >im 
moving  toward  the  flame  from  a  distance  of  0.5ft  with  a  velocity  of  14ft/sec.  Curves  are 
drawn  for  different  geometric  sizes  of  cylinder. 

4.  Results  And  Conclusions 

Knowing  B  the  droplet  vaporization  rate  can  be  found  by:  /}={&/ p  )(X/c  )//Kl+B  )  The 
results  are  shown  in  Fig  1  &  2.  An  interesting  feature  noticed  by  the  curves  in  Fig  2.  When 
the  droplet  temperature  after  reaching  a  maximum  value,  starts  to  decrease  as  droplet  gets 
closer  to  the  radiation  source.  This  is  not  noticed  in  the  pre-evaporation  zone  of  systems  of 
relatively  low  radiative  heat.  Because  droplet  heat  up  progresses  very  slowly  relatve  to  the 
residence  time  in  cold  zone.  If  a  very  long  cold  region  or  a  long  residence  time  is  taken,  the 
same  feature  namely  the  reduction  of  droplet  temperature  is  noticed.  The  authors  explanation 
for  this  is  that  total  radiation  heat  to  the  droplet  is  reducing  with  a  factor  proportional  to  d2. 
And  when  droplet  radius  is  very  small  radiation  to  droplet  gets  is  very  small  and  practically 
there  is  no  external  source  of  heat  to  supply  the  latent  heat  of  vaporization  and  droplet  uses 
its  own  thermal  heat  to  vaporize  and  hence  droplet  cooles  down  very  slowly.  In  models 
neglecting  radiation  heat,  and  only  considering  convection  to  the  droplets,  temperature  stays 
steady  after  reaching  to  a  maximum.  In  fact  as  convection  coefficient  h  is  inversly 
proportional  to  the  droplet  diameter(D).  Then  total  convective  heat  to  the  surface  of  the 
droplet  would  be  proportional  to  D.  But  total  radiation  to  the  droplet  is  proportional  to  D2, 
and  effect  of  droplet  diameter  is  more  significant  in  models  including  radiation  heat 

Fig.  3  shows  aseries  of  curves  for  radiation  flux  to  droplets  at  the  centerline  of 
cylindrical  chambers.  The  wall  radiation  becomes  important  in  many  cases.  For  example  in 
the  cylinders  of  30cm  radius  with  30cm  hot  zone  and  300cm  radius  with  150cm  hot  zone 
thickness,  the  wall  radiation  contributes  65%  and  45%  of  the  total  radiation  respectively.  This 
happens  when  the  wall  temperature  is  taken  to  be  the  same  as  hot  zone  temperature{2000  °C). 
It  is  interesting  to  point  out  that  the  variation  of  the  radiative  flux  in  space  upstream  of  the 
flame,  is  not  drastic.  As  a  first  approximation  constant  radiation  flux  can  be  considered  in 
some  engineering  calculations. 
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Abstract 

The  affects  of  different  gas  and  liquid  phase  models 
on  the  vaporization  behavior  of  a  single  component 
isolated  droplet  are  studied  for  both  the  stagnant  and 
the  convection  situations.  In  conjunction  with  four 
different  liquid  phase  models,  namely  d2-law,  infinite 
conductivity,  diffusion  limit  and  internal  vortex 
circulation,  the  different  gas-phase  models  include  a 
spherically  symmetric  model  in  the  stagnant  case  and 
ftanz-Marshall  correlation  plus  two  other  axisymmetric 
models  in  the  convective  case.  A  critical  comparison 
of  alt  of  these  models  is  made.  The  use  of  these 
models  in  a  spray  situation  is  also  examined.  A 
transient  one-dimensional  flow  of  an  air-fuel  droplet 
mixture  is  considered.  It  is  shown  that  the  fuel  vapor 
mass  fraction  can  be  very  sensitive  to  the  particular 
liquid-phase  and  gas-phase  models. 


Nomenclature 


■  specific  heat  at  constant  pressure.  cal/gm/ *K 

•  D  De,  gas  diffusivity,  cm2/sec 

•  H  C  '  T  \  heat  from  the  gas-phase  given 
to  tRe  liquid  phase 

per  unit  mass  of  fuel  vaporized,  cal/gm 

•  length  of  the  tube,  cm 

•  L  C  '  T  \  heat  of  vaporization,  cal/gm 

P  e 

■  L  tr  \  ratio  of  gas-phase  characteristic  length 
and°the  initial  drop  size 

■  Mf  Mc',molecuiar  weight  of  fuel 

•  Mq  Mfi*.  molecular  weight  of  oxidizer 

•  rJ’  vaporization  rate,  gm/sec 
*  r  c  o 

•  n^/L  2,  number  of  droplets  per  unit 
cross  section 

area.  1/cm2 

•  P  Pc'»  flea  pressure,  atmosphere 

■  rfc  r  \  droplet  radius,  n 

•  gas  temperatura  at  the  tuba  entranca.  *K 
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•  T  T  %  gas  temperature 

■  Tfc  T  *,  droplet  surface  temperature 

•  t  \c',  time 

»  t  /0  7L  2,  ratio  of  convection  time  to 
diffusion  time 
in  the  gas-phase 

■  gas  velocity  at  the  tube  entrance,  cm/sec 

•  v  vc  *  9®s  velocity,  cm/sec 

■  Vk  Vc',  droplet  velocity 

•  x  Le>  axial  coordinate,  cm 

■  Lc.  droplet  location 

•  fuel  vapor  mass  fraction  at  the  droplet 
surface 

■  fuel  vapor  mass  fraction 

■  oxidizer  mass  fraction 

■  neutral  species  mass  fraction 

•  spatial  grid  size 

■  ratio  of  specific  heats 

•  m  9®s  viscosity,  gm/cm/sec 

•  f  fc’,  variable  defined  in  Eq.  (11) 

■  p  pe\  gas  density 

■  liquid  fuel  density 

•  p  *lp  \  ratio  of  initial  gas-phase  density  and 
thce  liquid  fuel  density 


1.  Introduction 

Spray  vaporization  and  combustion  studies  are  of 
primary  importance  in  predicting  and  improving  the 
performance  of  systems  utilizing  spray  injection. 
Combustors,  fire  suppression,  spray  drying  and  various 
forms  of  chemical  power  plants  are  typical  examples  of 
these  systems.  Often,  the  vaporization  of  a  droplet  in 
the  spray  is  affected  by  neighboring  droplets. 
However,  in  spray  combustion  computations,  it  is 
assumed  that  the  overall  spray  behavior  cen  be 
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obtained  by  summing  behavior  of  individual  isolated 
droplets  surrounded  by  a  gas -phase  that  itself  has 
varying  properties.  Even  when  the  assumption  that 
droplets  behave  as  if  they  were  isolated  from  each 
other  is  not  satisfactory,  the  behavior  of  a  single 
isolated  droplet  in  an  oxidizing  environment  will  still 
provide  a  fundamental  input  to  the  overall  spray 
analysis. 

A  complete  kinematic  and  vaporization  study  of  an 
isolated  droplet  in  a  convective  field  involves  the 
simultaneous  description  of  the  transport  processes  in 
both  the  gas  and  the  liquid  phases,  the  coupling 
between  the  phases,  the  phase  change  processes  at  the 
interface,  and  the  chemical  reaction  in  the  gas-phase. 
However  such  a  detailed  model  is  too  complex  to  be 
used  in  the  spray  computations  at  this  time.  As  a 
result,  several  simplifying  assumptions  are  made  in  the 
spray  analysis.  For  example,  it  is  assumed  that  the 
gas-phase  processes  around  the  droplet  are  quasi¬ 
steady,  and  that  a  phase  equilibrium  exists  at  the 
interface  with  the  phase  change  being  given  by  a 
C  laus  i  us -Cl  ape  yr  on  type  relation.  Also  the  assumptions 
of  constant  gas-phase  transport  properties  and  of  a 
spherically  symmetric  flow  field  with  a  convective 
correction  are  often  employed.  In  addition,  the  liquid 
phase  transport  processes  may  be  completely  neglected 
by  using  a  d2-law  (constant  droplet  temperature)  model 
[1]  or  by  using  an  infinite  conductivity  (spatially 
uniform  liquid  temperature  in  the  droplet)  model  [2]. 
The  models  based  on  the  above  assumptions  are 
somewhat  useful  in  that  they  simplify  the  problem 
considerably  and  thus  provide  easier  solution  to  the 
problem.  However  there  is  a  need  to  upgrade  the 
droplet  models  in  the  spray  analysis  so  that  a  better 
physical  representation  is  given  to  transient  behavior, 
liquid  phase  heat  and  mass  transport,  internal 
circulation,  and  gas-phase  boundary  layers.  At  the 
same  time,  there  is  a  need  to  simplify  the  more 
complex  droplet  models  so  as  to  enhance  their 
usefulness  for  the  spray  analysis.  The  work  in  this 
paper  was  undertaken  to  fulfill  the  above  premise. 


The  besic  droplet  vaporization/combustion  model  for 
an  isolated  single  component  droplet  in  a  stagnant 
environment  was  given  by  Godsave  [1].  Since  then 
this  model  has  been  extensively  studied  both 
experimentally  and  theoretically.  These  studies  have 
been  reviewed  by  Williams  [4],  Faeth  [5]  and  Law  [8J. 
More  sophisticated  studies  have  also  been  reported. 
These  consider  the  effects  of  relaxing  the  restrictions 
of  the  basic  model.  Law  and  Law  [7]  considered  the 
variable  gas-phase  transport  property  effects.  Hubbard 
at  a)  [8]  considered  the  transient  effects  for  the 
spherically  *  symmetric  case  and  concluded  that  the 
quasi-steady  gas-phase  approximation  yields  small 
errors,  at  least  for  pressures  below  ten  atmospheres. 
Their  comparison  with  the  constant  transport  property 
models  using  reference  conditions  indicated  that  the  1/3 
rule  gives  the  best  agreement.  There  have  also  been 
studies  on  the  transient  droplet  heating.  These  include 
the  infinite  liquid  phase  conductivity  model  [2]  and 
conduction  limit  model  [9].  It  should  be  noted  that 
the  infinite  conductivity  model  is  an  artificial  one  for 
the  spherically  symmetric  case  since  no  internal 
mechanism  is  present  to  increase  the  liquid  heating  rate 
above  that  of  the  conduction  limit. 

In  many  practical  situations,  the  droplet  vaporizes  in 
a  convective  gas  field.  The  gas-phase  convection 
influences  the  vaporization  process  in  two  ways.  First 
it  increases  the  gasification  rate  as  well  as  the  heat 


transfer  rate  between  the  phases.  Secondly  it 
generates  liquid  circulation  inside  the  droplet  which 
increases  the  liquid  heat  transfer  rate.  Note  that  at 
very  high  relative  velocity,  droplet  deformation  and 
shattering  can  occur.  This  situation  is  not  considered 
in  this  paper.  Semi-empirical  correlations  [10]  account 
for  the  gas-phase  convection  by  expressing  the 
veporization  and  interphase  heat  transfer  rates  as  a 
modification  of  the  spherically  symmetric  case. 
Sirignano  [11]  analyzed  the  convective  case  through  a 
combination  of  stagnation  point  and  flat  plate  analyses 
and  concluded  that  the  convective  case  should  not  be 
treated  by  a  correction  on  the  spherically  symmetric 
case.  Prekash  and  Sirignano  [12]  analyzed  the  gas- 
phase  end  the  liquid-phase  flow  fields  for  a  single 
droplet  in  a  convective  gas  field.  They  considered  a 
gas  phase  boundary  layer  outside  the  droplet  end  a 
Hill's  vortex  in  the  droplet  core  with  thin  viecous  and 
thermal  boundary  layers  near  the  droplet  surface  end  an 
inviscid  internal  wake  near  the  axis  of  symmetry.  Their 
results  show  that  the  infinite  conductivity  case  is  never 
realized  and  the  characteristic  liquid  phase  heat 
diffusion  time  reduces  (from  the  diffusion-limit  case)  by 
an  order  of  magnitude.  Their  analysis,  however,  is  too 
complicated  and  computer-time  consuming  to  be 
included  in  spray  calculations.  Tong  end  Sirignano  [13] 
simplified  this  analysis  substantially  by  neglecting  the 
thin  boundary  layer  inside  the  droplet  surface.  The 
results  of  this  simplified  model  are  in  close  agreement 
with  those  obtained  from  the  more  exact  analysis  [12]. 
The  simplified  vortex  model  and  the  other  liquid  phase 
heating  models  for  a  single  droplet  and  for  e  spray 
situation  are  critically  axaminad  in  the  following 
sections.  In  Section  2,  the  effects  of  d2-law,  infinite 
conductivity,  conduction  limit,  end  vortex  models  on  the 
veporiratfon  behavior  of  a  single  component  droplet  ere 
studied.  Both  the  stagnant  and  the  convective 
environments  are  considered.  For  the  convective  case, 
Ranz-Marshall  correlation  [10],  an  optimum  average  of 
the  stagnation  point  region  and  the  flat  plate  (droplet 
shoulder)  region,  end  a  more  exact  axisymmetric  model 
ere  studied.  Afl  these  gas-phase  and  liquid  phase 
models  ere  examined  for  three  hydrocarbon  fuels, 
namely,  n-hexane.  n-decane.  end  n-hexedecane.  The  use 
of  different  liquid  phase  models  in  e  spray  vaporization 
situation  is  described  in  Section  3.  where  a  transient 
one-dimensional  flow  of  an  air-fuel  droplet  mixture  in 
an  open  tube  is  considered.  The  effects  of  different 
droplet  models  on  the  fuel  vapor  mass  fraction 
distribution  in  the  tube  ere  discussed  in  detail.  An 
unsteadiness  in  the  gas-phase  properties,  which  is 
entirely  due  to  the  discrete  droplet  group  locations,  is 
also  discussed.  Conclusions  ere  stated  in  Section  4. 

2.  Single  Droplet  Vaporization 

Basically,  the  existing  literature  on  single  droplet 

veporization  can  be  classified  into  two  major 

categories:  spherically  symmetric  end  axisymmetric. 

The  different  models  in  these  two  categories  ere 
discussed  below. 

2.1  Spherically  Symmetric  Models 

2.1.1  d2  Lsw  Model.  The  most  notable  earlier  work 
on  droplet  vaporization  is  by  Godtsve  [1].  in  thst 
study,  a  quasi-steady  sphericslly-symmetric  model  was 
used.  The  droplet  temperature  was  assumed  to  be 
uniform  and  remained  constant  at  its  wet  bulb  value. 
The  properties  in  both  the  gas  end  liquid  phases  were 
assumed  to  be  constant,  together  with  Lewis  number 

equal  to  unity.  At  the  gas-liquid  interface,  it  was 

eesumed  that  the  fuel  vapor  mats  fraction  was  e 
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function  of  the  surface  temperature  given  by  eome 
equilibrium  vapor  pressure  equation  such  as  the 
Cleusius-Clapeyron  relation.  This  theory  gives  the 
classic  d2Haw  and  is  the  simplest  possible  modal 
describing  droplet  vaporization.  It  should  be  noted  that 
this  model  neglects  the  liquid  phase  heat  and  mass 
transfer  and  is  basically  a  gas-phase  modal.  It  does 
not  consider  some  of  the  important  physics  and  only 
yields  a  crude  estimate  of  the  droplet  vaporization  rata. 

Some  of  the  resulting  relationships  ere  given  below; 


m 
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equations  (1)  and  (2)  remain  the  same  except  that  L  is 
replaced  by 


L 


m3 


The  additional  term  accounts  for  the  transient  liquid 
droplet  heating.  Subscript  l  denotes  the  liquid  phase. 

2.1.3  Conduction  Limit  Modal.  Whan  the  internal 
liquid  motion  is  not  significant,  heat  transfer  inside  the 
droplet  will  be  controlled  by  thermel  diffusion  only. 
This  will  bs  e  reasonable  model  for  the  stagnant  case 
and  represents  the  slowest  heat  transfer  limit 
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where  m  is  the  fuel  mass  vaporization  rata.  B  is  the 
Spalding  Transfer  number,  t  is  the  latent  heat  of  fuel 
and  is  the  fuel  vepor  mass  fraction  at  the  interface. 
Subscripts  g  and  s  denote  gas-phase  and  droplet 
surface,  respectively. 


Note  that  ■  y  (Tf)  it  a  function  of  surface 
temperature  only.  Therefore  T  and  B  can  be 
determined  by  equation  (2)  and  subsequently  can  be 
obtained  through  equation  (1).  For  a  given  ambient 
condition,  T#  and  B  are  fixed  in  this  model,  equation  ID 
can  be  written  ee 


—  •  k  •  constant  (3) 

dt  s 

which  says  that  the  radius  squared  (proportional  to 
tha  droplet  surface  areal  decreases  linearly  in  time  and 
hance  this  modal  is  raferrad  to  as  d2-law  model  in  the 
literature.  *  j 

2.1.2  Infinite  Conductivity  Model.  In  e  combustor,  the  [ 
droplet  is  initially  cold  and  heats  up  with  time.  The  ! 
initial  transient  droplet  heating  period  has  only  recently 
been  given  attention.  Lew  [2]  studied  droplet 
combustion  with  rapid  internal  mixing  where  droplet  1 
tempareture  is  spatially  uniform  but  varying  with  time. 

It  was  found  that  droplet  heating  is  a  significant  cause  1 
of  the  unsteadiness  of  droplet  combustion  and  should 
bs  taken  into  account  in  any  realistic  analysis  of j 
unsteady  droplet  combustion  phenomena.  Basically,  this  , 
model  is  tha  earns  as  the  d2Haw  model  except  that  the  j 
constant  liquid  phase  tempareture  assumption  is  relaxed  i 
and  la  replaced  by  s  uniform  but  time-varying  j 
tempareture  inside  the  droplet.  The  gas-phase  model 
remains  spherical ly  symmetric  end  quasi-steady. 


Some  authors  believed  that  this  uniform  temperature 
limit  is  related  to  the  rapid  internal  liquid  circulation 
limit  end  hence  is  referred  to  as  rapid  mixing  model. 
Sirignano  [11]  showed  that  even  in  the  limit  of  high 
vortex  strength,  the  internal  liquid  circulation  can  only 
reduce  the  character i at ic  langth  scale  for  diffusion  by 
an  order  of  magnitude.  The  rapid  mixing  limit  can 
never  exist.  Rather,  it  would  be  conceptually  more 
correct  to  think  that  the  uniform  temperature  limit 
results  from  the  infinite  conductivity  limit.  Hence  It  is 
more  appropriate  to  eell  it  an  infinite-conductivity 
modal. 


Tha  liquid  pheta  tamparatura  variation  is  governed  by 
the  well  known  heat  diffusion  equation. 


with  the  initial  end  boundary  conditions 


(I)  at  t  •  0;  T(r,  t)  ■  T  (r> 

O 


aT(r.t) 

(ii)  at  r  •  0;  — - —  •  0  (S) 

dr 

dT(r.t) 

(lit)  at  r  ■  r#(t);  k^  — - —  ■  q^(t) 


where  TQ(r)  is  the  initial  temperature  distribution  and 
q^(t)  is  the  liquid  phase  heat  flux  at  the  droplet  surface. 

Equations  (D  and  (2)  remain  essentially  the  same 
except  that  L  is  replaced  by 
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The  additional  term  accounts  for  the  liquid  phase  heat 
flux  at  the  surface. 


Since  the  droplet  is  vaporizing,  the  diffusion  fitld  has 
•  moving  boundary.  Obviously  this  model  is  more 
complicated  then  the  infinite  conductivity  limit  model. 

2J2  Droplet  Vaporization  with  Convective  Effects 

2^.1  Renz-Mershell  Model.  The  studies  discussed  so 
far  ere  baaed  on  the  model  of  spherical  symmetry,  in 
reality,  the  practical  problem  of  droplet  vaporization  in 
a  spray  involves  a  convective  situation  in  which  there 
is  a  relative  gas-dropiet  velocity.  Many  investigators 
[4.5]  suggested  empirical  correlations  for  vaporization 
rate  in  a  convective  field  as  s  correction  to  the 
spherical  symmetric  case.  The  typical  form  of 
correlation  is  m  •  m.  .  flRe,  Pr)  where  f(Re. 

Pr)  is  the  correction  rector.  ^Tn  tne  present  study,  the 
Ranz-Marshall  correlation  is  examined,  for  which  the 
factor  f(Re,  Pr)  is  given  by 

f(R»,  Pr)  -  1  ♦  0.3  R,1/J  Pr1'3  (6) 


Although  this  type  of  correlations  is  very  simple, 
there  is  really  very  little  theoretical  justification  for 


them.  Moreover,  tfiese  correlation*  ere  based  on 
experiments  conducted  under  quest-steady  conditions. 
As  mentioned  before,  there  is  a  transient  droplet 
heating  stage  and  during  that  period  the  quasi* 
staadiness  assumption  is  invalid.  Indeed  Sirignano  [11] 
showed,  through  theoretical  analysis,  that  a  correlation 
of  Renz-Marshall  type  cannot  satisfactorily  predict  the 
vaporization  rate. 

XX2  Prakash  and  Sirignano's  Axisymmetric  Model. 
Recently.  Prakash  and  Sirignano  [12,  14],  studied  the 
problem  of  transient  liquid  droplet  vaporization  in  a  hot 
convective  environment.  They  first  [14]  studied  liquid 
internal  circulation  and  droplet  heating  and  later  [12] 
developed  a  gas-phase  boundary  layer  analysis  end 
coupled  it  to  the  previous  liquid  phase  analysis. 
Basically,  they  used  a  two-dimensional  axisymmetric, 
model  and  divided  the  problem  Into  several  physical 
regions  as  shown  in  Figure  1.  These  regions  were: 
outer  inviscid  gas  flow,  gas-phase  boundary  layer,  liquid 
phase  boundary  layer,  internal  liquid  wake  and  inviscid 
liquid  core.  Because  of  the  relatively  high  Reynolds 
number  and  high  Peclet  number,  the  boundary  layer 
approximation  was  used  in  t heir  analysis  for  both  the 
momentum  and  the  energy  transfers  in  both  the  gas  and 
the  liquid  phases.  The  outer  inviscid  gas  flow  was 
treated  as  steady  potential  flow  around  a  sphere  with 
no  flow  separation. 

The  gas-phase  and  liquid-phase  boundary  layers  were 
both  treated  as  quasi-steady  for  both  momentum  and , 
energy  transfer,  while  for  the  liquid  core,  the 
momentum  transfer  was  treated  as  quasi-steady  and  the 
energy  transfer  was  considered  transient.  This  model.  j 
named  as  Prakash  and  Sirignano's  model  here,  has  been; 
included  in  the  single  droplet  study.  It  should  be 
mentioned  that  this  model  although  quite  detailed,  uses 
an  algorithm  that  is  too  cumbersome  to  be  included  in 
a  complete  combustion  analysis.  Simplifications,  taking 
into  account  the  important  physics,  are  needed. 

In  Prakash  and  Sirignano's  analysis,  there  is  a  thin 
thermal  boundary  layer  near  the  droplet  surface  which 
is  coupled  to  the  thermal  core  in  the  matching  region. 
The  thermei  boundary  layer  which  allows  for  the  heat 
flux  to  adjust  from  two-dimensional  behavior  along  the 
droplet  surface  to  one-dimensional  behavior  along  the 
boundary  layer  end  thermal  core  matching  region,  was 
treated  as  quasi-steady.  The  importance  of  the  thermal 
boundary  layer  and  the  quasi-steadiness  assumption 
have  recently  been  reviewed  by  Tong  end  Sirignano 
[13].  The  results  of  that  study  show  that,  unless  the 
thermal  boundary  layer  is  very  thin  (very  large  Peclet 
numbers),  the  thermei  inertia  term  is  important  and 
quasi-steadiness  assumption  for  that  layer  is  invalid. 
The  results  also  tend  to  suggest  that  the  elimination  of 
the  thermei  boundary  layer  representation  end  the 
assumption  that  the  thermei  core  solution  is  valid  up  to 
the  droplet  surface  may  still  give  solutions  with 
acceptable  degree  of  accuracy.  It  should  be  noted  that 
while  the  thermei  boundary  layer  formulation  is  two 
dimensional,  the  thermei  core  formulation  is  one-i 
dimensional.  Hence,  the  removal  of  the  thermal 
boundary  layer  wilt  simplify  the  problem  tremendously. 
With  the  elimination  of  the  thermei  boundary  layer  end 
the  assumption  that  ths  thermal  cor#  solution  being 
valid  up  to  the  droplat  aurlaca,  Tong  end  Sirignano 
[13]  subsequently  simplified  the  thermei  core  solution. 
Their  simplified  liquid  phase  model  is  included  in  this 
study. 


2.2.3  Tong  and  Sirignano's  Axisymmetric  Model.  Tong 
end  Sirignano  [15]  further  developed  a  one  dimensional 
gas-phase  model.  They  simplified  the  axisymmetric 
convective  analysis  in  the  gas-phase  by  representing 
the  heat  and  mass  transfer  rates  by  an  optimum 

average  for  the  stagnation  point  region  and  the  shoulder 
region  of  a  droplet.  The  results  from  coupling  this 

simplified  gas  model  to  their  simplified  liquid  phase 

model  which  wes  obtained  earlier  are  in  reasonably 
good  agreement  with  the  results  of  the  more  detailed 
modal  of  Prakash  and  Sirignano.  This  gas-phase  model 
has  been  employed  in  the  spray  calculations  to  be 
reported  in  Section  3. 

2.3  Results  and  Discussions 

Liquid  fuel  droplets  initially  at  300  K  vaporizing  in 

1000  K,  10  atm  pressura  fuel-free  environment  is  used 
to  study  the  various  gas-phase  and  liquid-phose  models. 
Hexane,  decane  and  haxadecane  droplets  of  initial  radius 
of  0.00476  cm  are  considered.  The  physical  properties 
are  the  same  as  those  in  the  Prakash  and  Sirignano's 
analysis  [12]. 

Tha  results  for  the  stagnant  case  ere  given  in  Figures 
2  end  3.  The  gas-phase  model  is  quasi-steady  and 
spherically  symmetric.  Three  different  liquid  phase 
models  namely  the  d2Haw,  the  conduction  limit  and  the 
infinite  conductivity  limit  are  compared.  Since  the 
existence  of  a  vortex  in  a  stagnant  environment  cannot 
be  raalizad  conceptually,  it  is  excluded  for  comparison. 
One  would  think  that  if  the  vortex  model  were  included, 
it  would  be  in  between  the  infinite  conductivity  limit 
end  the  conduction  limit.  This  indeed  can  be  shown  to 
be  true.  Figures  2  and  3  show  the  temporal  variation  in 
the  surface  temperature  and  surface  area  respectively. 
The  physical  time,  t  is  nondimensionelized  by  using  the 
liquid  thermei  diffusivity,  end  the  initial  droplet 
radius,  r^,  as  r  •  Tor  the  d2-iaw  model  the 

surface  temperature  is  at  the  wet-bulb  temperature  and 
remains  constant  while  the  surface  area  regresses 
linearly  in  time.  The  d2-lew  serves  as  an  asymptotic 
limit  for  tha  other  two  models.  This  is  reflected  in 
Figure  2  where  the  surface  temperature  of  the  other 
two  models  approaches  the  wet  bulb  temperature  and 
Figure  3  where  the  curves  become  more  linear. 

Since  the  d2-lew  neglects  the  initial  transient  droplet 
heating,  tha  droplat  vaporizes  much  fester  than  the 
other  two  models.  The  difference  is  bigger  for  the 
heavier  fuels  which  have  higher  boiling  temperatures 
and  hence  longer  transient  heating  period  Conversely, 
the  conduction  limit  end  the  infinite  conductivity  limit 
have  about  the  same  droplet  lifetime.  Note  that  these 
two  limits  intersect  each  other.  Since  the  droplet 
temperature  is  uniform  in  the  infinite  conductivity  limit, 
the  surface  temperature  increase  per  unit  of  energy 
absorption  is  lass.  Consaquantfy.  tha  fuai  vapor  mass 
fraction  at  tha  droplat  surfscs  is  lowar,  which  leads  to 
lower  vaporization  rate.  On  the  other  hand,  the 
difference  between  the  ambient  temperature  end  the 
surface  temperature,  which  relates  directly  to  the  heat 
transfer  rate  to  the  liquid,  is  higher.  This  results 
eventually  in  e  higher  surface  temperature  and  faster 
vaporization  rate  during  the  latter  portion  of  the 
lifetime  for  the  infinite  conductivity  model.  Note  that 
the  wet-bulb  temperatures  in  these  calculations  are 
considerably  lower  then  their  corresponding  boiling 
temperatures  of  the  fuels.  The  effect  of  the  different 
liquid-phase  models  on  the  interior  liquid  temperature  is 
quite  significant  and  therefore  the  conduction  limit 
model,  which  it  considered  to  be  exact,  should  be  used. 
If  detailed  temperature  distribution  inside  the  droplet  is 


needed.  The  infinite  conductivity  model  may  be  useful 
in  the  low  ambient  temperature  case  when  the  droplet 
lifetime  is  long.  The  d2-law  model  which  gives  poor 
agreement  with  the  conduction  limit  model  is 
oversimplified  and  should  be  discarded. 

The  results  for  the  convective  case  are  given  in 
Figures  4  to  7.  In  these  calculations,  the  Reynolds 
number  (based  on  droplet  radius)  is  initially  100  and 
decreases  with  the  radius  as  the  droplet  vaporizes.  The 
relative  velocity  is  2500  cm/sec  and  is  assumed  to  be 
constant  throughout.  In  reality  the  relative  velocity  will 
be  reduced  by  the  drag  force;  this  is  considered  later  in 
the  spray  vaporization  case. 

Although  Prakash  and  Sirignano's  model  is  the  most 
detailed  and  should  be  considered  to  be  the  most  exact, 
it  is  quite  cumbersome  to  be  included  in  a  complete 
spray  analysis.  Instead.  Tong  and  Sirignano's  model 
and  the  simplified  vortex  model  which  has  been  shown 
[15]  to  give  close  agreement  with  Prakash  and 
Sirignano's  results  are  used  as  reference  model  in  the 
figures. 

Figures  4  and  5  show  the  results  of  the  various  liquid 
phase  models.  The  d2Haw  does  not  exist  in  this 
convective  case.  Again,  the  conduction  limit  and  the 
infinite  conductivity  limit  intersect  due  to  the  same 
reasoning  as  that  in  the  stagnant  case.  The  vortex 
model  lies  in  between  the  other  two  models  as 
expected.  Note  that  the  droplet  lifetime  is  considerably 
shorter  than  that  in  the  stagnant  case.  In  the 
convective  case,  (r/r  )3/2  regresses  asymptotically 
linearly  with  time  and  it  is  used  in  Figures  5  and  7. 
Figure  5  shows  that  for  the  heavier  fuel  the  percentage 
variation  in  droplet  lifetime  can  depend  quite 
significantly  on  the  liquid  phase  model. 


later  at  the  tube  entrance  are  prescribed.  The  injection 
of  fuel  droplets  is  intermittent.  One  group  of  droplets 
is  injected  every  given  time  interval.  The  number  of 
droplets  in  a  group  (or  with  a  characteristic)  represents 
the  number  of  droplets  per  unit  cross-section  area.  The 
•  frequency  of  injection  and  initial  droplet  velocity 
depends  on  the  overall  fuel-air  ratio  in  the  tube,  the 
mass  flow  of  air.  droplet  spacing  and  droplet  diameter 
at  the  tube  entrance.  The  prescription  of  the  overall 
fuel-air  ratio,  the  droplet  diameter  and  the  droplet 
velocity  at  the  entrance  yields  the  value  of  the  droplet 
number  density.  The  droplet  spacing  in  the  axial 
direction  and  the  number  of  droplets  per  unit  cross- 
section  area  are  so  adjusted  as  to  provide  initially  an 
isotropic  droplet  spacing  in  a  unit  cube.  Initially  a  high 
relative  velocity  is  provided  between  the  two  phases. 
As  a  droplet  group  moves  in  the  hot  gas  stream,  it 
accelerates.  At  the  same  time,  the  droplets  heating  and 
vaporization  is  taking  place.  The  spray  processes 
influence  the  state  of  the  gas,  i.e.,  the  gas-phase  is 
continuously  retarded,  cooled,  and  enriched  with  fuel 
vapor.  The  gas-phase  properties  are  also  being 
influenced  by  the  upstream  convection.  All  these  gas- 
phase  and  liquid-phase  processes  are  modeled  by  a 
system  of  unsteady,  one-dimenstonal  equations.  The 
gas-phase  is  represented  in  Eulerian  coordinates 
whereas  the  liquid  phase  is  represented  in  Lagrangian 
coordinates.  The  governing  equations  in  the  non- 
dimensional  form  are  given  as  follows: 

GAS  PHASE  EQUATIONS 
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Figures  6  end  7  show  the  effects  of  the  different 
gas-phase  models.  The  Ranz-Marshall  correlation 
overpredicts  the  vaporization  rate  and  underestimates 
(for  decane)  the  surface  temperature.  As  indicated  by 
Sirignano  [11],  this  type  of  correlation  has  very  little 
theoretical  justification  and  cannot  give  satisfactory 
results.  The  present  study  supports  Sirignano's 
qualitative  analysis.  Tong  and  Sirignano’s  model  appear 
to  be  in  good  agreement  with  the  more  detailed  mode) 
of  Prakash  and  Sirignano. 

The  results  indicate  that  the  droplet  heating  and 
vaporization  are  essentially  unsteady  for  most  of  its 
lifetime,  particularly  for  the  heavier  fuel.  For  a  detailed 
analysis.  Prakash  and  Sirignano's  model  is 
recommended  but  for  more  practical  application,  Tong 
and  Sirignano's  model  with  the  simplified  vortex  model 
would  be  useful. 

X  SPRAY  VAPORIZATION  STUDY 


In  the  spray  vaporization  problem,  a  one-dimensional  < 
transient  flow  of  air  and  fuel  droplets  in  an  open  tube 
is  considered,  where  the  motion  and  vaporization  of  a, 
monodispersed  spray  in  a  laminar,  hot  gas  stream  are 
studied.  The  purpose  is  to  examine  the  effects  of 
different  liquid  phase  models  on  the  bulk  vaporization 
characteristics  as  well  as  on  the  gas  medium.  The 
different  liquid-phase  models  considered  are  the  infinite 
conductivity,  the  conduction  limit,  and  the  vortex 
models.  The  different  convective  models  for  the  gas- 
phase  are  the  Ranz-Marshall  correlation  and  the  model 
of  Tong  and  Sirignano.  The  physical  situation  consists 
of  a  continuous  laminar  flow  of  hot  air  in  an  open 
tuba.  Tha  gas-phase  properties  initially  in  the  tube  and 
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LIQUID-PHASE  EQUATIONS 
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and  the  reference  temperature  is  obtained  by  1/3  rule 
as  recommended  by  Sparrow  and  Gregg  [  18]. 
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for  Ram-Marshall  correction  model  (15) 


for  Tong  -  Sirignano  Model 
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The  drag  coefficient  is  evaluated  by  using  an 
expression  as  proposed  by  Putnam  [19].  The 
evaluation  of  Reynold's  number  for  the  drag  coefficient 
is  based  on  the  free  stream  density  and  viscosity 
evaluated  at  the  1/3  reference  state.  This  follows  the 
recommendation  of  Yuen  and  Chen  [20].  It  should  also 
be  noted  that  the  effect  of  relative  droplet-gas  velocity 
on  the  vaporization  rate  [see  Eq.  (15)]  has  been  treated 
by  a  semi-empirical  correlation  [10].  as  well  as  by  an 
axisymmetric  model  [15]. 

3.1  The  Solution  Procedure 


24  /  Re  2/3\ 


(17) 


2  r  p 

Re.  • - ^  I  V-V. 

k  trLf  ,  '  k 


(18) 


i  rb  r  .  3/2i 


(19) 


where  T 


raf 


1 

-  r 

3 


(20) 


The  source  terms  in  the  gas-phase  equations  result 
from  the  coupling  be  ween  the  two  phases.  These  are 
given  in  the  Appendix.  It  should  be  noted  that  a 
transformation  as  given  by  Eq.  (11)  has  been  employed. 
This  transformation  is  useful  for  a  constant  volume 
situation  [16].  However,  it  has  been  retained  here  (a 
constant  pressure  case)  for  the  sake  of  generality. 


In  non-dimensionaiizing  the  gas-phase  equations,  the 
length  scale  is  the  length  of  the  tube,  and  the  velocity 
scale  is  the  gas  velocity  at  the  tube  entrance.  The 
time  scale  is  determined  by  these  two  scales.  The 
gas-phase  properties  are  non-dimensionelized  by  using 
the  respective  properties  at  the  entrance.  For  the  liquid 
phase  equations,  the  droplet  location,  velocity  and 
surface  temperature  are  respectively  non- 
dimensionalized  by  the  gas-phase  length,  velocity  and 
temperature  scales.  The  droplet  radius  is  non- 

dimensionalized  by  the  initial  droplet  radius.  The  above 
non-dimensionalization  gives  rise  to  three  additional 
dimensionless  groups.  tr.  L,  and  pr;  tf  is  the  ratio  of 
convective  time  to  diffusion  time  in  the  gas-phase.  Lr 
is  the  ratio  of  gas-phase  length  scale  and  initial  drop 
radius,  and  pf  is  the  ratio  of  initial  gas-phase  density 
and  liquid  density. 


A  hybrid  Euierian  Lagrangian  numerical  scheme  is 
used  to  calculate  the  gas-phase  and  liquid-phase 
properties  at  the  (r»«T)m  time  level  from  the  values 
known  at  the  n*  time  level.  First  of  all.  a  second- 
order  accurate  scheme  is  employed  to  interpolate  the 
gas-phase  properties  from  the  Euierian  locations  (fixed 
grid  points)  to  the  Lagrangian  (or  droplet)  locations. 
The  scheme  uses  the  gas-phase  properties  at  two  grid 
points  x.  and  x.a^1  and  gives  the  corresponding 

I  properties  at  a  Lagrangian  location  X  ,  where  Xk  is 

I  between  x(  and  x.^.  Using  these  gas-pnase  properties, 
ithe  droplet  surface  temperature  is  then  calculated.  For 
ithe  infinite  conductivity  model,  it  involves  the  solution 
)  of  an  ordinary  differential  equation  for  each  droplet 

group.  For  the  conduction  limit  and  vortex  models,  a 
partial  differential  equation  needs  to  be  solved.  The 

1  details  have  been  discussed  in  Section  2.  Knowing  Tk 
)  for  each  droplet  g roup,  other  droplet  properties  (Xk,  V  . 
rk)  can  be  obtained  by  solving  Eqs.  (13)  -  (16).  A 
second-order  Runge-Kutta  scheme  is  used  for  this 
purpose.  It  should  be  noted,  however,  that  there  is 
some  loss  of  accuracy  due  to  the  fact  that  the  gas- 
phase  properties  are  not  updated  in  the  Runge-Kutta 

scheme.  Using  the  new  liquid-phase  properties,  the 

source  terms  as  given  in  the  Appendix  can  be  evaluated 
iat  the  Lagrangian  locations.  Then  a  second-order 
i  accurate  scheme  is  used  to  distribute  these  source 

I  terms  from  a  Lagrangian  location  to  the  two 

I  neighboring  gas-phase  locations.  Using  these  source 
j  terms,  the  gas-phase  properties  (Yf.  Yq  and  f)  at  the 
(n+D*  time  level  are  obtained  by  solvmg  Eq.  (8).  An 
explicit  finite-difference  method  is  used  for  this 
purpose.  The  gas  temperature  and  gas  density  are  then 
obtained  by  using  Eqs.  (11)  and  (12).  Finally  the  gas 
velocity  is  obtained  from  an  integral  form  of  the 
continuity  equation  (7).  It  is  worth  emphasizing  that  in 
the  present  case,  the  gas-phase  convection  term 
dominates  the  diffusion  term.  Therefore,  an  upwind 
difference  scheme  [21]  is  employed  for  the  convection 
term.  It  should  also  be  noted  that  an  explicit  method 
was  found  to  be  the  most  efficient  for  reacting  two- 
phase  flow  [21].  However,  an  implicit  method  may  be 
more  efficient  for  the  present  vaporizing  case  without 
combustion  since  the  equations  may  not  be  as  stiff. 


The  important  assumptions  made  in  writing  these 
equations  are  that  the  gas  pressure  is  constant, 
radiative  heat  transfer  is  negligible,  the  species 
diffusion  follows  Pick's  law  with  equal  mass 
diffusivities  for  each  pair,  the  specific  heat  at  constant 
pressure  is  constant,  and  the  gas  phase  Lewis  and 
Schmidt  numbers  are  unity.  In  addition,  the  product  p 
D(*  p)  is  assumed  constant.  It  is  noteworthy,  however, 
that  in  the  calculation  of  liquid-phasa  properties  p  is 
considered  a  function  of  temparature  as  given  by 
Sutherland  correlation  [17].  As  indicated  in  Eq.  (19),  a 
referenca  temperature  is  used  for  the  calculation  of  p 


3.2  Discussion  of  Results 

The  effects  of  using  different  liquid-phase  models  on 
the  bulk  spray  and  gas-phase  properties  are  now 
discussed.  The  various  values  used  in  the  computation 
are  listed  in  Table  1.  The  criterion  used  in  selecting 
these  values  was  to  consider  a  spray  vaporization 
situation  with  a  moderate  gas-phase  convection  for  the 
droplets.  The  selected  values  provide  an  initial 
Reynold's  number  (based  on  the  relative  velocity  and 
| the  drop  radius)  of  56.5  and  a  droplet  residence  time  of 
I  about  12J5  jnsec.  During  this  time,  eighty  percent  of 


the  droplet  mass  Is  vaporized.  The  effects  of  different 
liquid-phase  models  on  the  bulk  gas-phase  properties  is 
presented  in  Figs.  8-11.  Figures  8  and  9  indicate  that 
the  choice  of  a  particular  droplet  heating  model  can 
significantly  influence  the  fuel  vapor  distribution  in  the 
tube.  This  can  have  a  profound  effect  on  the 

subsequent  combustion  process  in  a  practical  situation. 
The  maximum  difference  between  the  fuel  vapor  mass 
fraction  values  is  as  much  as  one  hundred  percent 
between  the  conduction  limit  and  the  infinite 
conductivity  models*  As  Fig.  9  indicates,  the 
conduction  model  initially  (up  to  a  distance  of  5  cm) 
predicts  the  highest  fuel  vapor  mass  fraction,  whereas 
the  infinite  conductivity  model  predicts  the  lowest.  At 
later  times  the  situation  is  reversed.  The  vortex  modal 
results  are  always  in  between  those  of  the  other  two 
models.  The  difference  in  the  fuel  vapor  distribution  is 
a  direct  consequence  of  the  difference  in  the  droplet 
surface  temperature  values  for  the  three  models.  As 
discussed  in  Section  2,  the  conduction  limit  model 
initially  predicts  the  highest  droplet  surface  temperature 
and  therefore  the  fastest  vaporization  rata,  whereas  the 
infinite  conductivity  model  predicts  the  lowest  surface 
temperature  and  the  slowest  vaporization  rate.  As  also 
discussed  in  Section  2.  this  behavior  is  reversed  at 
later  times. 

The  fuel  vapor  mass  fraction  profiles  in  Figs.  8  and  9 
exhibit  an  oscillatory  behavior.  This  is  due  to  the 
periodic  nature  of  the  droplet  injection  process  which 
causes  a  finite  droplet  spacing  in  the  streamwise 
direction.  Thus  the  period  of  this  oscillation  is  the 
same  as  the  time  interval  between  two  subsequent 
injections.  This  is  clearly  illustrated  in  Fig.  10.  As 
this  figure  indicates  the  profiles  overlap  after  every  0.5 
msec,  which  is  the  time  interval  between  two 
subsequent  injections.  This  unsteady  gas-phase 
behavior,  which  is  entirety  due  to  the  spray 
discretization,  should  be  kept  in  mind  in  spray 
modeling.  The  gas  temperature  profiles  at  time  ■  16 
msec  for  the  three  liquid  models  are  shown  in  Fig.  11. 
The  continuous  drop  in  the  gas  tamperature  is  indicative 
of  the  continuous  cooling  of  the  gas-phase  due  to 
droplet  heating  and  vaporization.  The  difference  in  the 
gas  tamperature  distribution  due  to  the  different  liquid 
phase  models  is  not  as  significant  as  that  in  the  fuel 
vapor  distribution.  This  is  due  to  the  fact  that  the 
spray  haating  and  vaporization  constitute  only  a  small 
heat  loss  in  the  overall  gas-phase  heat  budget.  The 
oscillatory  behavior  in  tha  gas  temperature  profiles  is 
again  due  to  the  intermittent  droplet  injection  process. 
This  intermittency  is  prescribed  exactly  in  our  idealized 
calculation.  In  a  practical  situation,  it  may  not  be  so 
well-defined  and  perhaps  should  be  taken  as  random. 
That  is.  the  distance  between  droplets  as  they  flow 
into  the  domain  of  interest  will  not  be  fixed  as  they 
are  given  in  this  calculation  but  rather  will  be 
distributed  about  soma  average  value. 

The  effect  of  different  liquid  heating  models  on  tha 
droplet  properties  is  illustrated  in  Table  2.  The 

properties  of  the  first  droplet  group,  i.e.  the  location, 
tha  non-dimensional  radius,  the  surface  temperature,  and 
tha  velocity  are  given  as  functions  of  time  for  the 
three  models.  The  difference  in  the  surface 
temperature  and  droplet  radius  values  for  the  three 

models  is  very  similar  to  that  discussed  in  Section  2. 
The  conduction  limit  model  predicts  the  highest  droplet 
velocity,  whereas  tha  infinita  conductivity  model 

predicts  the  lowest  value.  This  occurs  since  the 
conduction  limit  model  initially  predicts  the  fastest 

vaporization  rata,  tha  smallast  drop  sizt.  and 


consequently  the  highest  droplet  acceleration.  It  should 
be  noted,  however,  that  the  difference  in  the  droplet 
velocities  for  the  three  models  is  quite  small. 

The  effect  of  different  convective  models  is 
presented  in  Figs.  12  -  13.  It  should  be  noted  that  the 
model  of  Tong  and  Singnano  is  valid  only  when  the 
Reynolds  number  is  significant;  as  indicated  by  Eq.  (15), 
this  model  predicts  a  vaporization  rate,  which  is 
proportional  to  the  square  root  of  the  Reynolds  number. 
For  this  reason,  initial  droplet  radius  of  100  microns  is 
used  for  these  calculations.  This  gives  a  value  of  112 
for  the  initial  Reynolds  number.  Figure  12  shows  that 
the  Ranz-Marshall  correlation  overpredicts  the 
vaporization  rate  as  compared  to  that  given  by  the 
axisymmetric  model.  This  behavior  is  consistent  with 
the  one  observed  for  the  single  droplet  case  in  the 
previous  section.  Figure  13  gives  the  differences  in  the 
gas  tamperature  distribution  for  the  two  models.  Since 
the  Ranz-Marshall  correlation  model  predicts  a  higher 
vaporization  rate  and  therefore  a  higher  heat  loss  from 
the  gas  phase,  it  gives  lower  gas  temperatures  than 
those  given  by  the  other  model.  It  is  also  important  to 
note  that  at  very  small  Reynolds  number,  the 
axisymmetric  model  predicts  an  unacceptably  small 
vaporization  rate.  Consequently  a  more  realistic  spray 
model  should  have  the  capability  to  switch  from  the 
axisymmetric  model  to  a  Ranz-Marshall  correlation  type 
of  model  as  the  Reynolds  number  becomes  small.  This 
is  presently  under  investigation. 


4.  Conclusions 

The  most  common  vaporization  models  have  been 
compared  for  fuels  of  varying  volatibility  and  droplets 
in  a  high-temperature  environment.  Both  isolated 
droplet  and  spray  vaporization  hava  been  studied  The 
use  of  the  d2-law  or  the  infinite-conductivity  model 
(sometimes  named  the  rapid  mixing  model)  has  been 
shown  to  be  very  inadequate.  For  sphencally- 
symmetric  vaporization  (no  relative  gas-droplet  motion), 
the  conduction  limit  model  for  transient  heating  is 
recommended.  In  practical  situations  where  a  relative 
gas-droplet  motion  exists,  the  simplified  vortex  model 
of  Tong  and  Sirignano  is  recommended  wfien  the 
Reynolds  number  based  on  relative  velocity  is  high. 
That  model  compares  well  with  the  more  detailed 
model  of  Prakash  and  Sirignano  and  predicts  well  the 
effects  of  the  laminar  gas-phase  boundary  layer  over 
the  droplet  and  the  internal  circulation  in  the  droplet. 

The  spray  calculations  indicate  an  inherent 
unsteadiness  due  to  the  intermittent  character  of  a 
spray.  This  unsteadiness  will  be  interesting  whenever 
we  wish  to  resolve  structure  on  the  scale  of  the 
average  distance  between  dropfets. 
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Appendix 

The  source/sink  terms  in  the  gas-phase  equations  (1. 
2)  are  given  as  follows: 
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In  the  above  relations.  Syf  and  Syo  are  tha  sourca  and 
sink  tarms  for  tha  fual  vapor  and  oxidizer  mass 
fractions,  xf^  is  tha  mola  fraction  at  tha  droplat 
surfaca,  M#'  is  tha  molecular  weight  of  tha  gas^phasa 
(excluding  fuel  vapor),  T^'  is  tha  liquid  fual  boiling 
temperature  at  tha  normal  pressure  p  p  is  tha 
prevailing  pressure,  and  j  represents  tfie  summation 
over  all  droplat  groups  which  happen  to  be  in  a  given 
gas-phase  mash,  of  size  Ax.  A  phase  equilibrium 
assumption  has  bean  made  at  tha  droplat  surface  and 
the  Clausius-Ctapeyron  relation  [Eq.  29]  has  been  used 
Tha  different  methods  of  computing  tha  droplat  surface 
temperature  T  distinguish  tha  various  liquid  phase 


models  examined  in  the  present  spray  problem, 
have  been  discussed  in  Section  2. 


These 


M6.J  FLOW  WGIOKS  OtfTS  IOC  MO  NlTHlN  A  FUEL  D**l£T  AT 
HIGH  REYNOLDS  NUMBER  (frm  Ref  22) 


net  'wti  ***o«  ruts  mactioh  *wwi  oi*i*»ct  *t 
«•  mu.i«eoi»*  ran  otrrintnT  uouio-phah  metis 


I 


i 


•  I  I'M*  iM 


rit.M  Ctl  ^kU  KfTfiATUW  »IMU*  OUTMICI  *! 

«  milimco**  ran  oirnnm  tJOuJD*i**Att  «oo*l* 


r 

r 

f 

!. 


a 

n 


h 
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Liquid  Density 
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Initial  Drop  Radius 
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Initial  Or  op  Temperature 
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NUMERICAL  MODELING  OF  ONE-DIMENSIONAL  ENCLOSED 
I  HOMOGENEOUS  AND  HETEROGENEOUS  DEFLAGRATIONS 

S.  IC  Aggarwal  and  W.  A.  Sirignano 
Mechanicai  Engineering  Department 
Carnegie-Mellon  University 
Pittsburgh,  Pennsylvania  15213 


ABSTRACT 

The  objective  of  this  study  is  to  evaluate  the  effectiveness  of  several  finite 
difference  techniques  in  modeling  unsteady  homogeneous  and  heterogeneous  flames 
in  closed  combustors.  The  problem  of  propagation  of  a  pre-mixed  homogeneous 
flame  and  a  heterogeneous  flame  in  air-fuel  droplet-fuel  vapor  mixture  is  formulated. 
The  governing  equations  have  been  solved  using  five  different  numerical  techniques, 
an  explicit  method,  an  iterative  implicit  method  without  any  formal  linearization,  an 
iterative  implicit  method  with  quasi-linearization,  a  non-iterative  time-linearization 
block  implicit  scheme,  and  an  operator-splitting  method.  The  results  are  obtained  in 
terms  of  the  profiles  of  gas  and  liquid  phase  properties  at  different  instances  of 
time  as  the  flame  traverses  the  combustor.  From  these  results,  the  relative 
efficiency  and  accuracy  of  each  of  the  numerical  techniques  are  discussed. 
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ABSTRACT 

A  hybrid  Eulerian-Lagrangian  method  is  employed  to  study  a  laminar  unsteady 
propagating  flame  through  an  air/fuel-vapor/fuel-droplet  mixture  in  a  one-dimensional 
combustor.  The  gas-phase  properties  are  obtained  by  using  a  continuum  approach, 
whereas  the  computation  of  liquid-phase  properties  is  based  on  a  Lagrangian 
approach.  A  quiescent  mixture  of  air  and  fuel  mist  is  ignited  by  providing  a  hot  wall 
at  the  end  of  the  combustor.  Following  ignition,  a  flame  initiates  and  propagates 
into  the  cold  spray  mixture.  The  characteristics  of  the  propagating  spray  flame  are 
presented  in  terms  of  the  profiles  of  gas-phase  and  liquid-phase  properties.  Effects 
of  droplet  heating  models,  overall  equivalence  ratio,  fuel  volatility  and  initial  droplet 
size  are  discussed.  Results  for  the  infinite-conductivity  mode)  and  the  conduction 
limit  mode!,  which  represent  the  fastest  and  the  slowest  heat  transport  modes  in  the 
droplet  interior  respectively,  indicate  that  the  droplet  heating  models  have  only  a 
minor  effect  on  the  bulk  spray  flame  characteristics,  when  the  overall  mixture  is 
close  to  stoichiometric.  However,  when  the  overall  mixture  is  either  fuel  rich  or  fuel 
lean,  the  differences  in  heating  models  are  significant.  An  increase  in  fuel  volality 
and/or  a  reduction  in  initial  drop  size  enhances  the  speed  of  flame  propagation,  and 
imparts  a  premixed  type  character  to  the  flame. 


